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SECTION  I 
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INTRODUCTION 

The  main  goal  of  this  study  of  High  Power  Laser  Coupling  was  to  examine 
the  coupling  of  laser  energy  to  solid  targets.  To  accomplish  this  goal, 
a broad  spectrum  of  experiments  was  carried  out  at  2.8  and  3.8  micron 

wavelengths  using  the  Boeing  pulsed  chemical  laser  ^and  at  10.6  microns 

(2' 

using  the  Marx  Bank  CO2  laser.  ' During  the  course  of  this  study,  the 
energy  output  of  the  chemical  laser  was  increased  from  152  to  292  joules  at 
HF  wavelengths.  In  addition,  operation  at  DF  wavelength  was  developed  yielding 
an  energy  of  144  joules.  By  operating  the  laser  at  various  combinations  of 
gas  pressure  and  flashlamp  voltages,  laser  effects  have  been  studied  at 
pulse  lengths  from  3 to  14  microseconds. 

The  laser  described  in  Section  II  was  used  in  an  extensive  series  of 
plasma  threshold  experiments  described  in  Section  III.  Seven  materials 
were  studied  and  thresholds  were  determined  for  HF  and  DF  wavelengths.  The 
materials  include  Lucite,  Grafoil , aluminum,  titanium,  stainless  steel,  nickel, 
and  copper.  In  general,  it  was  found  that  the  longer  wavelength  DF  radiation 
yielded  threshold  significantly  lower  than  HF  radiation.  Theoretical 
descriptions  of  plasma  production  thresholds  were  developed  in  terms  of 
simultaneous  fluence  and  irradlance  requirements  as  well  as  simultaneous 
bulk  temperature  and  irradiance  requirements. 

Section  IV  describes  the  results  of  thermal  coupling  experiments  at 
both  HF  and  DF  wavelengths  for  several  materials.  In  addition,  a limited 
number  of  experiments  were  carried  out  over  a range  of  ambient  pressures. 

In  almost  all  cases,  it  was  found  that  the  gross  thermal  coupling  efficiency 


1.  Hall,  R.  B.,  Maher,  W.  E.,  McClure,  J.  D. , Nichols,  D.  B.,  and  Pond,  C.  R. , 
"Laser  Beam  Target  Interaction,"  Vol . I,  2.8  Microns  Wavelength,  AFWL- 
TR-25-342,  July  1976. 


2.  Maher,  W.  E.,  Hall,  R.  B.,  and  Johnson,  R.  R.,  "Experimental  Study  of 
Ignition  and  Propagation  of  Laser-Supported  Detonation  Waves,"  J.  Appl. 
Ph/s.  45,  2138,  1974. 
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followed  the  well-known  trend.'  * ' The  coupling  coefficient  is  nearly 
constant  for  low  irradiances  and  then,  at  the  plasma  threshold,  jumps  in 
value.  Above  threshold  there  is  a broad  range  of  energies  for  which  the 
total  absorbed  energy  at  the  target  is  constant. 

Section  V describes  impulse  production  and  mass  loss  for  both  metal  and 
ablative  targets.  It  was  found,  as  expected,  that  lucite  (like  other  materials) 
has  a threshold  below  which  zero  mass  is  lost.  The  maximum  specific  mass  was 
found  to  be  0.18  gm/kJ  and  then  decreases  slowly  as  the  incident  fluence 
increases.  An  interesting  aspect  of  impulse  production  was  found  in  this 
series  of  experiments  when  an  I/E  value  of  75  dyne-sec/ joule  was  produced 
with  a transparent  target  backed  by  an  absorbing  layer.  Since  most  I/E 
ratios  are  less  than  10,  one  notes  that  large  impulse  production  occurs  any 
time  a nearly  confined  plasma  is  produced. 

Section  VI  discusses  new  work  on  thermal  coupling  coefficients  at  C02 
wavelengths.  By  using  both  different  sizes  of  targets  and  an  external  air 
flow,  the  plasma  coupling,  due  to  plasma  spreading,  was  determined.  The 
width  of  this  spreading  was  determined  as  a function  of  the  incident  laser 
energy. 

Section  VII  is  a tneoretical  study  of  melting  and  vaporization  of  targets 
of  aluminum,  titanium,  and  molybdenum.  Curves  of  melt  thickness,  melting 
rate,  and  impulse  production  were  obtained  as  a function  of  absorbed  lasnr 
irradiance.  It  is  found  that  high  melting  efficiencies  can  be  obtained  if 
a small  amount  of  vaporization  is  allowed.  Under  this  condition,  the  vapor 
pressure  produces  a radial  outward  force  which  removes  the  molten  material. 

It  is  found,  however,  that  melt  removal  by  pulsed  laser  beams  appears  to  be  only 
marginally  useful. 


3.  Rudder,  R.  R.  and  Augustoni,  A.  L.,  "Thermal  Deposition  Experiments  with 
Microsecond  Duration  CO?  Laser  Radiation,"  Laser  Digest,  AFWL-TR-75-229, 

1975. 

4.  Hall,  R.  B.,  Wei,  P.  S.  P.,  arid  Maher,  W.  E.,  "Laser  Beam  Target  Inter- 
action," Vol.  II,  Laser  Effects  at  10.6  Microns,  AFWL-TR-75-342,  July 

1976. 
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SECTION  II 


PULSED  CHEMICAL  LASER  CHARACTERISTICS 


1 . BACKGROUND 

In  order  to  carry  out  laser  effects  studies  over  a wide  range  of  para- 
meters, the  performance  of  the  Boeing  pulsed  chemical  laser  was  improved. 

This  laser  (PH0CL-10)  was  previously  described;^  a schematic  of  the  device 
is  shown  In  figure  1.  It  includes  a mixing  assembly  for  hydrogen  a.:d  fluorine, 
and  is  initiated  by  U-shaped  flashlamps  which  are  internal  to  the  reactor 
volume.  Buffer  cells  isolate  the  cavity  optics  from  the  reacting  gases. 

Both  the  flashlamp  current  pulse  and  the  laser  output  pulse  are  typically 
3 to  4 psec  full  width  at  half  maximum  (FWHM).  The  laser  energy  h*.o  reached 
152  J,  corresponding  to  3.5%  electrical  efficiency  before  starting  the  present 
study. 

2.  REACTOR  MODIFICATION 

A uv-reflecti ve  teflon  wall  lining  was  installed  in  the  laser  reactor 
to  increase  the  efficiency  of  fluorine  dissociation  by  flashlamp  radiation 
in  the  wavelength  band  235  t.o  335  mm.  Figures  2 to  4 show  the  variation  of 
several  performance  characteristics  as  a function  of  input  flashlamp  energy, 
using  the  teflon  covered  wall  with  only  two  flashlamps.  Figure  2 shows  a 
maximum  output  energy  from  two  lamps  equal  to  that  previously  obtained  with 
eight  lamps  with  no  uv  reflector.  The  gas  mixture  is  not  identical,  but  a 
Targe  part  of  this  improvement  from  3.5%  to  14%  electrical  efficiency  is  due 
to  the  Improved  uv  utilization. 

Two  gas  conditions  are  shown  in  the  figures.  The  circles  correspond  to 
a 1 atm  mixture  containing  10%  fluorine;  the  squares  are  for  a 0.6  atm 
mixture  with  7.4%  fluorine.  The  hydrogen  fraction  was  1.2%  in  either  case. 

The  lower-pressure  mixture  was  motivated  by  an  Interest  In  controlling 
pulse  length.  Pulse  lengths  are  displayed  In  figure  3.  It  is  clear  the 
conditions  giving  longer  laser  output  pulses  also  give  reduced  output  energy. 
This  result  Is  plotted  explicitly  in  figure  4 where  one  obtains,  for  example, 
a pulse  width  of  5 usee  with  80  J (joules)  of  laser  energy. 
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Reactor j Buffer 

Figure  1.  Diagram  of  Photolnitlated  HF-DF  Pulsed  Chemical  Laser 


ENERGY 

(JOULES) 


FLASHLAMP  ENERGY  (JOULES) 


Figure  2.  Variation  of  Laser  Output  as  a Function  of 
Flashlamp  Erargy  In.  Circles  Correspond 
to  1 atm  Initial  Pressure  Containing  10% 
Fluorine.  Squares  are  for  0.6  atm  With 
7.4%  Fluorine.  Hydrogen  Fraction  Is  1.2% 
In  Either  Case. 


Table  1 shows  several  benchmark  performance  points  obtained  during  the 
laser  development  phase  of  this  work.  The  highest  energy  point  of  figure  2 
Is  entered  as  the  first  line  of  table  1.  Chemical  and  electrical  efficiencies 
are  denoted  by  nc  and  ne  respectively. 

The  dependence  of  pulse  length  on  the  pressure  of  the  gas  mixture  Is 
shown  In  figure  5 over  the  range  0.2  to  1.0  atm  for  the  mixture  of  line  1, 
table  1.  Figure  6 shows  that  over  most  of  this  range,  0.3  to  1.0  atm,  the 
pulse  width  varies  approximately  as  the  Inverse  pressure. 

3.  EFFECT  OF  DILUENT 

The  Introduction  of  argon  as  diluent  resultea  In  a considerable  extension 
of  several  laser  performance  characteristics.  Line  2 of  table  1 shows  a 
factor-of- two- Improvement  In  electrical  efficiency  for  argon  compared  with 
line  1 for  nitrogen  diluent.  The  associated  differences  In  electrical 
properties  of  the  two  gases  Imposed  some  restrictions  on  the  laser  operating 
conditions  available,  but  argon  generally  allowed  the  wider  range  of 
performance.  The  improved  efficiency  can  be  attributed  In  part  to  the  lower 
heat  capacity  of  argon  and  resulting  differences  In  temperature-dependent 
rates  In  the  H2-F2  reaction  chain.  An  optimized  electrical  efficiency  of  33* 
Is  shown  In  line  3 of  table  1,  with  an  output  energy  over  100  J. 

The  effort  to  control  the  laser  puise  length  with  nitrogen  diluent  was 
continued  with  argon  diluent.  A condensed  version  of  the  results  of  this 
study  Is  shown  In  figure  7,  which  Includes  both  nitrogen  and  argon  diluent 
results.  Comparison  of  the  general  trend  of  argon  data  (triangles)  with 
nitrogen  data  (circles)  shows  that  for  a given  pulse  length  larger  energy 
per  pulse  can  be  obtained  with  argon.  Within  each  category,  the  solid  symbols 
of  figure  7 show  the  result  of  variation  In  flashlamp  voltage  and,  therefore, 
flashlamp  energy  and  level  of  Initial  photodissociation  of  fluorine.  The 
open  symbols  show  results  of  varying  the  Initial  gas  mixture  pressure.  The 
nitrogen  data  (open  circles)  with  Increasing  pulse  width  represent  total 
pressures  decreasing  from  1.0  to  0.2  atm.  The  argon  data  (open  triangles) 
with  Increasing  pulse  length  represent  pressures  of  1.0,  0.8,  and  0.6  atm. 

The  longest  pulse  operation  at  14  usee  and  70  0 Is  shown  as  line  4 of  table  1. 
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TABLE  1 


SELECTED  PERFORMANCE  CHARACTERISTICS  OF  THE  PHOTOINITIATED  PULSED  CHEMICAL 
LASER  WITH  THE  TEFLON -COVERED  WALL.  THE  HOUR  FRACTION  OF  H2  OR  D*  IS  1 .2% 
AND  THE  MIXING  TEMPERATURE  IS  AMBIENT.  THE  TWO  OUTPUT  COUPLERS  ARE  SAPPHIRE 
(87*)  AND  SILICON  (54%). 


1 Index  Lasant 

1 No- 

Diluent  p(total) 
(atm) 

°2 

F2 

No.  of 
Lamps 

Mirror 

Coupling 

E(laser) 

(J) 

FWHM 

(nsec) 

n n 

e c 
(*)  (*) 

1 

HF 

N2 

1.0 

0.25 

10 

2 

87 

148 

2.9 

14 

3.3 

2 

HF 

At 

1.0 

0.6 

10 

2 

54 

149 

3.5 

29 

3.3 

3 

HF 

Ar 

1.0 

o.e 

iO 

2 

54 

109 

5.0 

33 

2 1 

4 

HF 

Ar 

0.6 

0.5 

7. 

4 2 

54 

72 

13.8 

21 

2.7 

5 

DF 

Ar 

1.0 

0.6 

10 

4 

87 

144 

3.8 

14 

4.0 

» 

HF 

Ar 

1.0 

0.7 

10 

4 

87 

292 

2.7 

29 

8.1 

IS 


5 

TOTAL  f-RESSURk  (ATM) 


Figure  5.  Pulse  Length  Variation  with  Initial  Sas  Pressure 
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LASER  ENERGY  (JOULES) 

Figure  7.  FWHM  Pulse  Widths  Available  Using  Nitrogen  or  Argon  as  Diluent 


4.  OF  AND  HF  ENERGIES 

Efforts  to  introduce  a OF  output  capability  In  this  device  began  with 
nitrogen  diluent  with  disappointing  results.  The  substitution  of  argon  as 
diluent,  however,  gave  an  order  of  magnitude  Improvement.  Further  measure- 
ments led  to  the  conditions  listed  in  line  5 of  table  1,  with  144  J energy, 
using  four  flashlamps  and  the  teflon  wall  reflector. 

For  HF  wavelengths,  the  use  of  the  teflon  wall  reflector,  argon  diluent, 
and  four  flashlamps  led  to  the  optimum  energy  conditicns  listed  in  line  6 of 
table  1.  The  electrical  efficiency  is  29%,  and  the  chemical  efficiency  is 
8%.  For  both  HF  and  DF  wavelengths,  this  laser  produces  the  largest  energy 
per  pulse,  the  highest  electrical  efficiency,  and  the  highest  chemical 
efficiency  of  any  photolnitiated  chemical  laser  reported.  The  device  and  its 
development  are  reviewed  in  reference  5;  the  spectral  content  of  the  HF 
beam  is  discussed  in  reference  6. 

5.  BEAM  CHARACTERISTICS 

The  far-field  HF  beam  profile  was  measured  for  the  four- 1 amp  laser 
configuration  using  the  teflon  wall  lining  and  argon  diluent.  The  measure- 
ments were  made  in  the  focal  plane  of  the  f = 2 ni  mirror. 

The  general  shape  of  the  focused  beam  spot,  traced  from  the  burn  spot 
recorded  on  exposed  Polaroid  film.  Is  shown  as  the  dark  curve  in  figure  8a. 
This  pattern  is  a contour  of  equal  fluence  of  about  3 j/cm  , or  about  10* 
of  the  maximum.  The  four  Indentations  in  the  pattern  correspond  to  the  four 
flashlamps  In  the  reactor  and  their  insulating  feed-throughs.  An  analytic 
approximation  to  the  burn  pattern  Is  shown  as  the  light  curve  in  figure  8a. 

In  polar  coordinates,  this  shape  is  given  by 


5.  Nichols,  D.  B. , Hall,  R.  B.,  and  McClure,  J.  D.,  "Photoinitiation  F2  + H2/D9 
Chain-Reaction  Laser  with  Hlqh  Electrical  Efficiency,"  J.  Appl.  Phys.  47, 
4026,  1976. 

6.  Pond,  C.  R. , Hall,  R.  B.,  and  Nichols,  D.  B.,  “HF  Laser  Spec*  a!  Analysis 
Using  Near-Field  Holography,"  Appl.  Opt.  16,  67,  1977. 
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(1) 


r(8)  = ap  - b cos  (26) 


where  a„  = 1.73  cm  and  b = 0.57  cm. 
o 

Detailed  beam  scans  were  carried  out  in  this  same  focal  plane  using  a 
calorimeter  apertured  to  0.25  cm.  Figure  8b  shows  schematically  the  three 
directions  of  scanning,  illustrating  the  approximate  size  of  the  aperture 
and  spacing  of  the  data  points  on  the  beam  spot.  The  results  of  these 
measurements  are  shown  in  figure  9.  The  high-intensity  part  of  the  beam 
profile  is  flat-topped  as  found  previously  for  an  eight-lamp  configuration. 
The  two  scans  which  intersect  the  lamp  notches  show  a reduced  beam  width  as 
expected  from  the  burn  pattern.  The  largest  FWHM  (for  the  45°  scan)  is 
2.8  cm,  giving  a far-field  beam  divergence  of  14  milliradians. 

A "top-hat"  approximation  to  the  beam  spot  may  be  obtained  either  from 
the  actual  burn  shape  of  figure  1 or  from  equation  1.  The  radius  rQ 
of  this  equivalent  circle  is  related  to  the  45°  scan  radius  aQ  by  the  factor 


fe  = ro/ao 


where  fQ  < 1 due  to  the  areas  removed  by  the  lamp  images.  Direct  graphical 
analysis  of  the  burn  area  gives  f0  * 0.89.  From  equation  1,  we  obtain 


V2 


IT 

sT  l 


aQ  - b cos  4»  | d<J> 


which  gives  a similar  value,  fg  » 0.92. 

As  noted  above,  and  from  figure  8,  the  radius  of  the  focal  spot  Is  mainly 
characterized  by  Its  45°  parameters.  Thus,  as  noted  from  equation  2,  the 
H0-averaged"  radius  rQ  Is  approximately  0.89  times  the  45°  radius.  We  now 
define  the  equivalent  "top-hat"  radius  r0  by 


1— fj 


(4) 


*r8I  emax  = 2"  I r e{r)  dr 
0 

where  e(r)  is  the  energy  density  profile  at  45°  shown  in  figure  9 and  emax 

is  e(0).  The  integral  in  equation  4 was  evaluated  and  rfi  - 1.38  cm  obtained. 

The  overall  effective  radius  of  the  laser  beam  r,  is  then  defined  by 

2 

rl  = ^0re  suc^  that  the  incident  laser  energy  = tt  r^  emax- 

Table  2 summarizes  the  above  defined  beam  spot  characteristics  for 
three  mirror  focal  lengths  used  in  this  report. 


TABLE  2 

SPATIAL  PROFILE  PARAMETERS  OF  THE  FOCUSED  HF  LASER  BEAM  FOR  MIRRORS 
OF  VARIOUS  FOCAL  LENGTHS  F.  THE  VARIABLE  6 REFERS  TO  EQUATION  1. 

ALL  VALUES  ARE  IN  CM. 


Beam  Spot  Radius 

f = 200 

f * 76.2 

"*l 

II 

CO 

• 

<Sl 

P.adius  a^ 

(0.1  max,  0 * 45°) 

1.8 

0.69 

0.34 

Radius  rQ 

(r- Integrated,  0 = 45°) 

1.38 

0.53 

0.26 

Radius  r^ 

(r- Integra ted,  0-integrated) 

1.23 

0.47 

0.23 

HF  and  DF  laser  pulse  time  profiles  are  shown  in  figure  10  and  11, 
respectively.  Spikes  are  not  a prominent  feature  of  either  profile.  The 
HF  pulse  Is  quite  symmetric  with  typical  FWM  of  3.7  usee.  The  DF  profile 
has  a similar  rise  time  but  slightly  longer  duration,  with  FWHM  typically 
4.7  usee.  These  conditions  generally  apply  to  the  chemical  laser  results 
In  this  report  except  for  studies  which  explicitly  Incorporate  laser  pulse 
length  as  a variable. 
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SECTION  III 


THRESHOLD  CONDITIONS  FOR  FORMATION  OF 
TARGET  PLASMAS  BY  HF  AND  DF  LASER  RADIATION 


1 . INTRODUCTION 

With  increasing  exposure  of  a target  material  to  laser  radiation,  a 
point  is  reached,  in  terms  of  both  accumulated  energy  deposition  or 
instantaneous  Intensity,  at  which  a plasma  is  formed  at  the  target  surface. 

Under  these  conditions,  the  plasma  can  become  opaque  to  the  incident  radiation, 
shifting  the  principle  site  of  energy  deposition  from  the  target  surface 
to  the  plasma  Itself.  Because  this,  in  turn,  sustains  the  plasma,  one  forms 
"laser-supported  absorption  (LSA)  waves."  when  the  absorption  length  of 
the  plasma  Is  short,  and  the  plasma  front  propagates  up  the  incident  laser 
beam  at  supersonic  speeds,  these  waves  are  called  "laser-supported  detonation 
(LSD)  waves." 

Specification  of  the  conditions  for  the  formation  of  surface  plasmas 
is  of  considerable  interest,  due  to  the  wide  variety  of  effects  associated 
with  these  events.  These  effects  Include  the  enhancement  of  thermal  coupling 
of  the  incident  beam  into  the  target  material ,3^  increased  impulse  imparted 
to  the  target, and  the  Ionization  of  both  target  and  atmospheric  species. ^ 

Although  the  formation  of  surface  plasma  is  accompanied  by  several 
observables  In  both  the  target  and  surrounding  media,  definitions  of  the 
threshold  condition  have  been  largely  confined  to  two  categories.  In  the 
first,  visible  light  from  the  plasma  Is  observed,  often  by  fast  image 


7.  Ralzer,  Y.  P. , "Heating  of  a Gas  by  a Powerful  Light  Pulse,"  Soviet 
Physics  JETP  21,  1009,  1965. 

8.  Lowder,  S.  E.  and  Pettlnglll,  L.  C.,  "Double-Pulse  Pressure  and  Impulse 
Measurements,"  Lincoln  Laboratory  Optics  Research,  No.  2.,  1973,  p.  29, 
April  24,  1974. 

9.  Plrrl,  A.  N.,  "Theory  for  Momentum  Transfer  to  a Surface  with  a High- 
Power  Laser,"  Phys.  Fluids  16.  1435,  1973. 


90 


(2) 

conversion.'  ' In  the  second,  the  developing  opacity  of  the  plasma  to  the 
Incident  radiation  Itself  Is  detected.  A third  approach  has  been  to 
detect  plasma  formation  by  Its  electrical  conductivity.^^ 

2.  THEORETICAL  THRESHOLD  CONDITIONS  FOR  IGNITION  OF  LASER-SUPPORTED 
ABSORPTION  WAVES 

In  this  section  we  discuss  the  conditions  necessary  for  the  production 
and  maintenance  of  dense  plasmas  near  target  surfaces.  It  is  assumed  that 
two  simultaneous  conditions  must  be  satisfied  in  order  to  ignite  an  LSA 
wave.<"> 

a.  Simultaneous  Irradlance  and  Fluence  Requirements  (e+q  Model) 

First,  it  is  assumed  that  the  ignition  of  LSA  waves  occurs  by  heating 
of  defects  that  are  found  in  the  surface  of  typical  targets.  In  this  case, 
the  threshold  condition  is 


ChVet  joules/cm2  (5) 

where  h is  the  thickness  of  the  thin  defect,  C is  the  heat  capacity  of  the 
target  material,  and  Tv  is  the  vapor  temperature  corresponding  to  the 
ambient  pressure. 

It  Is  similarly  assumed  that  the  laser  irradlance  at  the  surface  qQ 
must  be  greater  than  a threshold  value  qt  in  order  to  maintain  an  LSA  wave. 

q0  > qt  (6) 


Let  us  first  assume  the  laser  pulse  Is  a step  function  of  time  of 
length  tp  as  shown  In  figure  12.  For  a focal  spot  area  S and  incident 


10.  Walters,  C.  T.,  Barnes,  R.  H.,  and  Beverly,  R.  W. , "An  Investigation  of 
Mechanisms  of  Initiation  of  Laser-Supported  Absorption  (LSA)  Waves," 
Battalia  Columbus  Laboratories,  Final  Report  on  Contract  DAAH01-73-A- 
0776,  1975. 

11.  Thomas,  P.  D.  and  Musal,  H.  M. , "A  Theoretical  Stu<(y  of  Laser-Target 
Interaction,11  LKSC-D352890  Lockheed  Missiles  and  Space  Group  Report, 
August  1973. 
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Figure  12.  Ratio  of  Fluence  at  Ignition  Tire  to  Incident  Fluence  as  a Function  of  E./Se 


laser  energy  E^,  the  laser  intensity  is 


% 


(7 


We  define  the  ratio  of  total  fluence  incident  on  the  target  at  the  time 
of  LSA  wave  Ignition  to  the  total  delivered  fluence  as 

fi  = (8) 

The  additional  requirement  that  the  laser  Intensity  be  greater  than 
the  threshold  Irradiance  leads  to 

Ei  > s vP  (9) 

If  q.t„  < e*  the  curve  of  as  a function  of  E,  is  shown  as  the 
t p t i 1 

g < 1 curve  of  figure  12,  where  the  parameter  g is  defined  as 

qtt 

9 * (10) 

et 

For  g < 1,  the  threshold  fluence  curves  versus  energy  are  again  shown  in 
figure  12.  In  general  (for  E » E^/Se^ 

fi  - 1 . 0 < E < 1 

- 1 , 1 < E < g 

» 1/E,  E < g (if  g <\) 


For  laser  pulse  shapes  other  than  square,  the  situation  is  more  complex. 
As  an  example,  we  next  consider  a laser  pulse  shape  given  by  (see  figure  13) 


for  0 < t < 2t„.  The  peak  laser  irradlance  q„  occurs  at  t = t„. 

p o p 


qo 


(13) 


The  irradiance  and  fluence  threshold  requirements  can  be  written 


(3/4)  E (2x-x2)  > g 

(14) 

(3/4)  E (x2-x3/4)  > 1 

(15) 

where  x * t/tp. 

When  both  of  the  conditions  specified  by  equations  14  and  15  are  met,  the 
ratio  of  absorbed  fluence  to  Incident  fluence  (Se/E^)  is  given  by 


- fj  - 3/4  (x2-x3/3)  (16) 


For  all  cases,  the  threshold  requirement  for  small  x is  determined  by  equation 
15.  For  larger  x,  equation  14  controls  the  threshold  requirements.  The 
boundary  xQ  between  these  two  regions  Is 

xo  ■ 7 [ 5 + 9 ' V’  ] (,7) 
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Figure  13.  Ratio  of  Fluence  at  Ignition  Time  to  Incident  Fluence  as  a Function  of  E./Se 


Thus,  we  find  the  conditions 


E = 


0 < x < xQ 
0/17  * (f°r  x > 1) 


3(x  -x  /3) 


(18) 


and  no  ISA  waves  are  possible  for  x > xQ  (when  xQ  > 1).  When  xQ  c 1 , one 
finds 


c - 4 0 < x < x„ 

c ~ o o ~ » 0 

3(x  -x/3) 


4q  x_  < x < 1 

Or  ’ 0 


(19) 


3(2x-x  ) 


and  no  plasma  waves  are  ignited  for  x > 1 (when  xQ  < 1). 

The  procedure  for  determining  curves  of  f,.  vs.  EfSe^  is  based  on  first 
choosing  a value  of  g and  then  computing  xQ  from  equation  17. 


(1)  For  xQ  > 1 


and  where 


^et  1 

0 < « f0 

. 3 2 1 3> 

(x~  " I xo  } 


ro  ' ? 


) 


(20) 


and  » 1 otherwise. 


(2)  For  xQ  < 1 


*i  • '/£*  0 < f1  < f( 


For  x < x < 1,  one  obtains  from  equation  14 


x ■ 1 


-fir 


(21) 
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For  a given  value  of  E^/Se^.,  a value  of  x is  computed  and  from  this 
value  f^  is  obtained  from  equation  16. 


It  is  easily  seen  from  equation  17  that  for  g < 1.5  that  xQ  > 1 . In 
this  case 


0 < E < 


T 


= 1/E 


E > 


(22) 


For  g > 1.5,  xq  < 1 and  one  finds 


f1  = 1/E  , 


(23) 


finally,  for  g > 1 .5 

f1  = 1 , E < y g (24) 

The  above  results  are  plotted  in  figure  13. 

b.  Simultaneous  Bulk  Temperature  and  Irradiance  Requirements  (T+q  Model) 

If  the  target  surface  Is  smooth  (free  of  defects)  or  if  the  wave- 
length is  short,  one  may  have  a threshold  condition  which  depends  on  the 
surface  temperature  itself.  Under  the  assumptions  of:  (1)  a laser  pulse  shape 
given  by  equation  12,  <2)  a surface  heating  rate  proportional  to  incident 
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irradlance  q,  arid  (3)  one-dimensional  heat  diffusion  Into  an  Infinitely  thick 

(121 

target,  the  surface  temperature  as  determined  from  Carslaw  and  Jaeger'  ' 

Is 


T(t)  „ = VkC^  x3/2  0“2x/5)  Agrees  K (25) 

where  x = t/t  , a is  the  material  absorptivity  at  the  laser  wavelength,  K is 
the  target  thermal  conductivity,  and  C is  the  target  heat  capacity.  For 
t > 2tp  (x  > 2),  the  temperature  declines  continuously.  The  maximum  surface 
temperature  as  found  from  equation  25  always  occurs  at  x = 3/2  . 

It  is  now  assumed  that  simultaneous  temperature  and  irradiance  thres- 
holds q^  must  be  surpassed  in  order  to  ignite  an  absorbing  LSA  wave.  These 
conditions  may  be  stated 


(26) 


In  equations  26  and  27,  Q is  the  ratio  of  the  laser  peak  irradiance  qQ  divided 
by  a target- dependent  irradiance  qy. 


12.  Carslaw,  H.  S.  and  Jaeger,  J.  C.,  Conduction  of  Heat  in  Solids,  Clarendon, 
Oxford,  1959. 
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Q " q0/qv 


h = V% 


qv 


I 


(28) 


It  Is  possible  to  consider  qy  as  an  equivalent  laser  irradiance  threshold 

required  to  raise  the  target  to  the  threshold  temperature  in  time  tp. 

For  an  aluminum  target  with  a = .06  and  t = 3.5  microseconds,  one  obtains 
7 2 P 

qy  = 3.4  x 10  watts/cm  . 

The  threshold  conditions  given  by  equations  26  and  27  are  analagous 
to  the  relationships  indicated  previously  by  equations  14  and  15.  It  is 
evident  from  the  laser  and  temperature  time  variation  that  threshold  require- 
ments are  not  met  for  x > 3/2.  It  should  also  be  noted  that 


Q - 


which  Is  then  of  the  form  given  previously.  The  dependence  of  Q on  x 
follows  equation  26  for  small  x and  then  for  x > xQ  follows  curve  27. 
The  boundary  x is  determined  by  equating  the  two  expressions. 


o -f  V 


(30) 


Table  3 shows  the  dependence  of  xQ  on  h obtained  from  equation  30. 


TABLE  3 


VARIATION  OF  BOUNDARY  x0(xQ-t/tp)  BETWEEN 
INTENSITY  AND  TEMPERATURE  CONTROLLED  THRESHOL  c 


.2 

4.37 

.3 

3.53 

.4 

3.01 

.6 

2.38 

.9 

1.81 

1.1 

1.53 

1.3 

1.28 

1.5 

1.021 

Curves  of  the  laser  energy  threshold  fi  as  a function  of  Q are  obtained 
in  different  regions  of  Q as  shown  in  Figure  14. 

For  h < 1.02  and  Q > 1.36,  the  curve  of  ff  vs.  Q is  obtained  by  the 
simultaneous  solution  of  equations  26  and  16.  This  threshold  curve  is  then 
dependent  only  on  target  vaporization  of  bulk  material.  For  Q < 1.36, 
fi  » 1.0. 


For  1.36  < Q < 1.67,  a discontinuity  in  fi  occurs  at  Q * h at  which  the 
value  of  fi  drops'- from  1.0  tc  the  above-determined  curve  for  h < 1.02. 

For  Q > 1.67,  fi  - 1 for  0 < Q < h.  At  Q = h,  a discontinuity  occurs 
reducing  fj  to  C.5.  For  h < Q < Q0  the  curve  of  fi  vs.  Q is  determined  by 
simultaneous  solution  of  27  and  16.  (The  value  Q0  is  that  value  of  Q 
occurring  at  x * xQ.)  For  Q > Qo»  fi  again  follows  the  h < 1.02  curve. 

3.  EXPERIMENTAL  METHOD  * . 

In  the  present  work,  plasma  formation  threshold  Was  determined  by 
observing  the  cutoff  of  Incident  laser  radiation  transmitted  tthrough  a 
small  hole  In  the  target. ^ The  main  elements  of  the  experimental  arrange- 
ment are  shown  schematically  in  figure  15.  Two  sample  beams  of  about  3X 
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Figure  14.  Ratio  of  Fluence  at  Ignition  Time  to  Incident  Fluence  as 
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Figure  15.  Diagram  of  Target  Room  Optics  and 
Experimental  Arrangement 


each  are  taken  from  the  primary  chemical  lase.'  beam  which  Is  incident  from 
the  left  In  figure  15.  One  of  these  beams  Is  directed  to  a calorimeter 
(Scientech  Model  3604).  The  other  beam  Is  focussed  by  a 165  cm  focal 
length  mirror  Into  a reflective  diffuser  and  a Ge(Au)  detector,  giving  a 
laser  pulse  time  profile  representative  of  the  entire  beam  cross  section. 

The  main  beam  goes  on  to  the  target  focussing  mirror. 

A second  calorimeter  was  placed  at  the  target  position  and  was  used  to 
calibrate  the  energy  monitor.  Beam  paths  from  the  CaF2  beam  splitter  to 
the  monitor  calorimeter  and  to  the  target  position  are  set  approximately 
equal  so  that  small  dally  variations  in  the  atmospheric  water  vapor  con- 
tent and  hence  In  the  atmospheric  absorption  coefficient  for  HF  laser 
radiation  are  compensated.  DF  absorption  over  these  distances  Is  negligible. 
The  optional  beam  path  to  the  vacuum  sphere  Is  also  shown  and  was  used  for 
the  work  reported  in  Section  V. 

For  the  present  threshold  measurements  targets  were  situated  at  the 
focus  of  a 76-cm  focal  length  mirror  and  were  oriented  at  8°  to  the  Incident 
beam  normal  to  prevent  feedback  from  specularly  reflective  targets  into  ti.e 
laser  and  calorimeter.  For  large  f -number  mirrors  this  effect  can  also 
distort  the  far-fleld  spatial  profile  by  redistributing  the  energy  at  the 
mirror  focal  plane.  The  beam  spatial  profile  Is  described  in  Section  II, 
and  has  an  effective  diameter  of  0.94  cm  for  this  mirror.  The  HF  and  DF 
far-fleld  profiles  are  nearly  Identical.  The  HF  time  profile  shown  in 
Section  II  has  3.7  psec  FWHM;  the  DF  pulse  Is  about  25<  longer. 

The  seven  materials  tested  included  five  metals  plus  Grafoll  and 
Luclte.  The  metal  samples  were  discs  with  2.22  cm  diameter  and  0.13 
cm.  thickness.  The  Grafoll  samples  were  0.025  cm  thick  and  2.54  cm 
square.  The  Luclte  samples  were  2.54  cm  squares  with  0.33  cm  thickness. 

All  materials  were  tested  In  the  "as  received"  condition;  the  only  surface 
preparation  was  cleaning  with  methanol  before  exposure.  A new  target 
was  used  for  each  laser  shot. 
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The  target  arrangement  for  threshold  measurements  Is  shown  schematically 
in  figure  15.  Pinholes  0.7  mm  In  diameter  were  centered  In  each  target. 

A short  aluminum  light  pipe  led  from  this  hole  to  an  attenuating  diffuser 
and  a second  Ge(Au)  detector  with  characteristics  matching  those  of  the 
other  fast  detector.  During  the  exposure  of  the  target  to  a laser  pulse, 
intensity  time  profiles  were  simultaneously  recorded  as  illustrated  In 
figure  16,  for  (1)  the  full  beam  (lower  trace)  and  (2)  the  laser  radiation 
transmitted  through  the  pinhole  (upper  trace),  which  is  truncated  by  the 
formation  of  an  LSA  wave. 

4.  THRESHOLD  DEFINITION  AND  ALUMINUM  RESULTS 

A threshold  time  ti  is  defined  as  the  time  at  which  the  transmitted 
intensity  falls  to  50%  of  the  uninterrupted  value,  and  the  fraction  of 
the  total  pulse  energy  (or  fluence)  delivered  before  t^  is  designated  fi; 
fi  is  referred  to  as  the  “transmitted  fraction"  of  the  incident  beam 
fluence 

Thresholds  will  be  discussed  in  terms  of  three  quantities:  (1)  total 
pulse  incident  energy  Ei  is  the  total  energy  directed  toward  the  target. 

(2)  Target  fluence  et  is  the  time- integrated  fluence  reaching  the  target 
before  cutoff,  or  E^fi/S  where  S is  the  beam  spot  area.  (3)  Maximum  target 
irradlance  c^is  the  largest  irradiance  which  reaches  the  target.  This 
maximum  occurs  at  ignition  time  ti  if  ti  is  before  the  pulse  peak  time  tQ; 
if  ti  > tQ»  q*  is  simply  the  peak  irradiance  of  the  incident  pulse. 

We  have  shown  previously  ^ that  threshold  values  of  each  of  these 
three  quantities  depend  strongly  on  whether  LSA  ignition  occurs  early  or 
late  In  the  pulse.  Thus  each  of  these  thresholds  Is  a function  of  ti,  or 
more  conveniently  f^.  Thresholds  previously  reported  with  no  reference  to 
an  associated  fi  are  ambiguous. 

Figure  16  shows  the  two  detector  time  profiles  for  a 44-J  DF  laser 
pulse  on  A1  6061;  the  transmitted  fraction  ^ ■ 0.71.  In  figure  17  the 
laser  pulse  energy  Is  Increased  to  83  J,  advancing  the  plasma  ignition  to 
an  earlier  time,  giving  fi  ■ 0.15.  Figure  18  Is  a plot  of  fi  versus 


Figure  16.  Incident  (Lower)  and  Transmitted  (Upper)  OF  Laser  Intensity  Profiles  for  E.  = 44  J 
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Figure  17.  Incl 


Incident  energy  for  both  HF  and  DF  radiation  on  aluminum.  Smooth  curves  are 
added  to  indicate  the  general  trends  of  the  data  without  regard  to  particular 
threshold  models. 

Two  definitions  of  threshold  energy  can  be  readily  derived  from  this 
graph.  The  first  is  the  value  associated  with  f^  * 1.0,  which  is  the 
largest  energy  deliverable  to  the  target  without  any  plasma  ignition  (as 
observed  by  cutoff  of  the  pinhole  transmission).  With  a finite  number  of 
points  available,  nearby  data  may  be  extrapolated  to  fj  * 1.  A second 
threshold  energy  Is  that  for  which  fi  = 0.5,  when  Ignition  occurs  halfway 
through  the  Incident  laser  pulse. 

Table  4 compares  HF  and  OF  thresholds  for  f^  = 1 , f i * 0.5,  and  f^  = 0.2 
for  A1  6061.  Also  listed  is  the  fluence  reaching  the  target  et»  again  as 
a function  of  f-j . The  maximum  intensity  reaching  the  target  % is  tabulated 
for  the  same  conditions.  It  appears  from  table  4 that  neither  e*  or  qm 
is  an  fj-independent  threshold  parameter  adequately  describing  all  the 
observed  thresholds  for  a given  material,  wavelength  band,  and  beam  spot 
size.  The  maximum  DF  irradiance  reaching  the  target  qm  is  found  to  be 
approximately  a constant  for  aluminum  but  not  for  other  materials. 

TABLE  4 

THRESHOLD  CONDITIONS  IN  TERMS  OF  INCIDENT  ENERGY  Ei , FLUENCE  REACHING 
TARGET  et,  AND  MAXIMUM  TARGET  IRRADIANCE  qm,  FOR  THREE  VALUES  OF  THE 
FRACTION  fi  OF  INCIDENT  HF  OR  OF  LASER  RADIATION  TRANSMITTED  THROUGH  A 
PINHOLE  IN  AN  ALUMINUM  6061  TARGET,  FOR  A 0,94  CM  DIAMETER  BEAM  SPOT. 


Lasant 

fi  « 1.0 

fi  = 0.5 

fj  S!  ( 

Ei  In  J 

HF 

45 

S3 

107 

DF 

49 

49 

65 

2 

et  In  J/cm 

HF 

€5 

60 

31 

DF 

?1 

35 

19 

qp,  In  MW/cm2 

HF 

18 

32 

33 

DF 

15 

15 

18 
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The  data  of  figure  18  are  replotted  In  figure  19  as  fluence  reaching 
the  target  e*  versus  Incident  fluence  e^.  The  HF  curve  shows  that  the 
maximum  et  deliverable  to  an  aluminum  target  (69  J/cm  ) occurs  for  fj  * 0.75 
and  does  not  appear  explicitly  in  table  4.  The  lines  through  the  origin 
indicate  loci  of  values  for  which  f^  * 1.0  and  f^  s 0.5. 

Figure  20  shows  the  aluminum  data  of  figure  18  compared  to  fits  using 
the  simultaneous  threshold  moo?^  described  in  Section  I II. 2.  The  e + q 
Model  curves  are  shown  as  solid  lines  while  the  T + q Model  calculations 
are  dotted  lines.  Equation  12  was  fit  to  the  time  profiles  of  figures  10 
and  11  to  give  tp  (HF)  = 2.6  gsec  and  tp  (DF)  « 3.3  ysec. 

For  DF,  the  curves  are  identical  above  fi  a 0.29  and  either  model  can 
provide  an  acceptable  fit  to  the  data.  Neither  model,  however,  adequately 
fits  all  the  HF  data.  The  HF  curves  shown  were  fit  to  the  higher  f^ 
points.  The  divergence  between  model  and  measurements  increases  with 
increasing  Ei  and  is  more  pronounced  for  the  e + q Model  curve. 


Similarly,  e + q Model  fits  have  been  generated  for  the  other  target 
materials  disciissedahjipvk  Most  of  these  fits  are  less  significant  than 
those  for  a 1 Jf  c sr  ,cewer  data  points  are  available;  these  fits 

are  not  display|w..  :ipf 

5.  THRESHOLD^^f  tiitNjjrfON  TARGET  MATERIAL 

Figures  21  througjB Ihow  the  transmitted  fraction  fi  as  a function 
of  Incldent^ncjjf/,  ircZrsurTi  for  both  HF  and  DF  wavelengths  for  six  other 
materials.  These  data  distributions  (and  the  associated  smooth  curves) 
were  Intended  particularly  to  give  good  approximations  to  the  f^  = 1 and 
fi  - 1/2  threshold  conditions.  Carbon  Is  included  as  Grafoil,  an  exfolia 
laminar  form  of  graphite.  Lucite  was  Included  to  allow  comparison  of 
the  metal  results  with  an  ablative  material . 
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Figure  19.  Fluence  Reaching  an  Aluminum  Target  Versus  Incident  Beam  Fluence 
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Figure  20.  Predictions  of  the  e+q  Model  (Solid  Lines)  and  T+q  Model  (Dotted  Lines)  Compared  to  the 
Aluminum  Threshold  Data  of  Figure  18.  Circles  are  HF,  Triangles  DF  Data. 


Figure  22.  Fluence  Fraction  Transmitted  for  Stainless  Steel. 
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Figure  26.  Fluence  Fraction  Transmitted  for  Luclte.  Circles  are  HF,  Triangles  DF  Data 


The  resulting  threshold  values  of  Incident  energy  , fluence  reaching 
the  target  et.  and  maximum  Intensity  are  tabulated  in  table  5. 

The  aluminum  results  of  table  4 are  Included  for  completeness. 

Estimated  uncerta Ini  ties  (standard  error)  in  the  threshold  values  are 
also  Indicated  In  table  3.  Errors  due  to  systematic  uncertainties  are 
estimated  as  15*.  This  value  governs  the  total  assigned  uncertainty 
except  In  those  cases  for  which  the  distribution  of  data  contributes  an 
additional  error;  the  values  Ei  (ft)  for  f]  * 1 and  1/2  for  Grafoll,  for 
example  (figure  26) , are  uncertain  by  + 20*;  this  is  combined  with  the 
systematic  error  to  give  the  +25*  Indicated  in  table  5. 

The  threshold  fluence  results  of  table  5 are  shown  In  graphical  form 
in  figures  27  and  28.  The  results  for  maximum  irradlance  are  shown  in 
figures  29  and  30. 


TABLE  5 

THRESHOLD  CONDITIONS  FOR  FORMATION  OF  LASER-SUPPORTED  ABSORPTION  WAVES  AS 
DETERMINED  BY  TARGET  PINHOLE  TRANSMISSION  CUTOFF,  FOR  SEVEN  MATERIALS 


f1  ■ 

Total 

Incident 

Pulse 

Enurgv 

E1 

(J) 

1 1/2 

Total 

Fluence 

Reaching 

Target 

et 

(J/cm2) 

1 1/2 

Maximum 

Irradlance 

qm  „ 
(MW/cnr ) 

1 1/2 

Assigned 

Uncertainty 

(*) 

Aluminum  HF 

45 

83 

65 

60 

18 

32 

25 

DF 

49 

49 

71 

35 

15 

15 

15 

Titanium  HF 

35 

58 

51 

42 

14 

23 

25 

DF 

30 

51 

43 

37 

9 

16 

25 

Stainless  HF 

45 

67 

65 

48 

18 

26 

15 

Steel  DF 

35 

56 

51 

41 

11 

17 

25 

Nickel  HF 

65 

82 

94 

59 

25 

32 

15 

DF 

30 

65 

43 

47 

9 

20 

25 

Copper  HF 

76 

83 

no 

60  ' 

30 

32 

15 

DF 

32 

36 

46 

26 

10 

11 

15 

Grafoll  HF 

65 

105 

94 

76 

25 

41 

20 

DF 

30 

66 

43 

48 

9 

20 

25 

Lucite  HF 

69 

175 

100 

127 

27 

69 

20 

DF 

69 

■*  — 

100 

-- 

21 

-- 

15 

69 


FLUENCE  REACHING  TARGET  (J/crn-*) 

Figure  27.  Comparison  of  Threshold  Fluences  for  f.  = 1 For  HF  and  DF  Laser  Radiation 


MAXIMUM  IRRADIANCE  (MW/cm*) 

Figure  29.  Comparison  of  Maximum  Irradiances  for  f.  = 1 for  HF  and  DF  Laser  Radiation 


SECTION  IV 


THERMAL  COUPLING  OF  HF  LASER  RADIATION  TO  METAL  TARGETS 

In  order  to  examine  the  efficiency  with  which  laser  energy  is  transferred 

to  target  materials,  a series  of  experiments  was  carried  out  at  high  Intensity 

in  ambient-pressure  air.  This  topic  is  of  Importance  because  of  the  well- 

known  enhanced  coupling  effect  In  which  plasma  Is  formed  above  a target.  The 

resulting  plasma  absorbs  the  laser  energy  and  then  radiates  optical  or  uv 
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radiation  toward  the  target.  ' For  reflective  targets  such  as  aluminum, 
the  effective  coupling  coefficient  can  be  increased  by  a factor  of  ten. (13,14) 


1.  THERMAL  COUPLING  COEFFICIENTS 

The  total  thermal  coupling  coefficient  a was  measured  for  four  metals 
as  a function  of  incident  HF  laser  energy.  The  coupling  was  measured  for 
two  beam  spot  sizes,  and  studies  were  made  of  the  effect  of  laser  pulse 
length  and  ambient  air  pressure  at  the  target. 

a.  Method 

Thermal  coupling  coefficients  of  the  metal  targets  were  obtained  by 
direct  calorimetry.  Chromel-alumel  thermocouples  were  formed  by  spot 
welding  the  0.013  cm  diameter  wires  to  the  target  rear  surfaces.  Axial 
diffusion  times  were  only  a few  milliseconds  due  to  the  nominal  0.13  cm 
target  thickness.  Radial  equilibration  times  were  determined  by  simultan- 
eous observation  of  time-dependent  temperature  readings  from  thermocouples 
mounted  at  the  center  and  at  the  edge  of  the  target.  The  time-dependent 
cooling  curves  from  center-mounted  thermocouples  were  then  extrapolated 
to  zero  time  to  obtain  an  effective  whole-target  temperature  rise  aT. 

The  associated  target  heat  was  divided  by  the  Incident  beam  pulse  energy 
to  give  the  thermal  coupling  coefficient  a. 


13.  Marcus,  S.,  Lowder,  J.E. , Manllef,  S.K. , and  Mooney,  D.L.,  'Laser 
Heating  of  Metallic  Surfaces,"  Project  Report  LIP- 31,  Lincoln 
Laboratory , May  1976. 

14.  McKay,  0.,  "High  Intensity  Coupling  Phenomena,"  Third  DOD  High  Power 
Laser  Conference,  Colorado  Springs,  November  1976. 
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b.  Thermal  Coupling  at  Large  Fluences 

The  variation  of  thermal  coupling  with  incident  laser  energy  was  studied 
over  a wide  range  of  fluence  (0  - 1300  J/cm2)  by  using  a 37.5  cm  collecting 
mirror  to  give  a beam  spot  of  effective  diameter  0.46  cm  as  discussed  in 
Section  II.  Target  diameters  for  this  work  were  1.59  cm.  Since  target 
heating  due  to  plasma  spreading  can  ocCur  well  outside  the  beam  spot, 
selected  high-energy  data  poii  ts  were  duplicated  with  larger  targets  (2.22  cm) 
to  show  that  the  target  heating  was  not  limited  by  target  size  for  these 
conditions. 


The  coupling  is  '■•letted  in  figure  31  for  aluminum  6061  as  a function  of 

incident  Hr  laser  energy.  Incident  pulse  fluence  and  peak  intensity  are 

^ 2 

also  shown.  The  ir.'adiance  scale  is  linear  up  to  about  220  MW/cm  due  to 

a constant  pulse  length  of  4.2  psec  FWHM.  Above  140  J energy,  modified 
gas  mixtures  gave  pulse  lengths  decreasing  to  3 usee  at  220  J. 


The  low-energy  data  of  figure  31  show  about  6%  absorbance  as  expected 
while  the  onset  of  plasma  formation  gives  a factor-of-three  enhancement,  in 
thermal  coupling.  The  half-height  point  in  this  enhancement  step  occurs  at 
about  120  J/cm^  fluence.  In  comparing  this  threshold  with  the  larger  spot  j * 
pinhole  thresholds  of  Section  III,  it  is  not  clear  what  value  of  f^  should  £ 
be  used.  However,  this  120  J/cm2  agrees  well  with  the  aluminum  results 
in  table  4 which  show  an  incident  fluence  of  2 x 60  J/cm  for  fj  = 0.5. 


The  solid  curve  of  figure  31  is  the  equation  a = (6.1  J)/E^.  This 

1 p 

simple  curve  describes  the  data  very  well  over  the  range  200  to  1400  J/cm  . 
The  result  is  a constant  absolute  energy  deposition  in  the  sample,  6.1  J, 
independent  of  incident  energy. 


Figures  32  to  34  give  similar  data  for  titanium,  stainless  stee^,  and 

nickel.  The  general  features  of  the  aluminum  results  are  repeated,  with  a 

2 

constant  energy  deposition  above  200  J/cm  . 
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Figure  3h  Thermal  Coupling  of  HF  Laser  Radiation  on  Aluminum 
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Figure  32.  Thermal  Coupling  of  HF  Laser  Radiation  on  Titanium 
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Figure  33.  Thermal  Coupling  of  HF  Laser  Radiation  on  Stainless  Steel 


68 


0 


200 


L 

o 


100 

ENERGY  ON  TARGET  (JOULES) 

I 1 » » « ‘ 1 « ■ » « I 

0.5  1.0  1-5 

FLUENCE  OCJ/CM2} 


-I I I I I I I 

100  200  300  400  600 

PEAK  IRRADIANCE  (MW/CM2) 


Figure  34.  Thermal  Coupling  of  HF  Laser  Radiation  on  Nickel 
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The  thresholds  observed  by  enhancement  of  thermal  coupling  are  compared 
with  pinhole  thresholds  In  table  6.  The  choice  of  * 0.5  type  pinhole 
thresholds  from  table  5 gives  good  agreement  with  the  thermal  enhancement 
thresholds. 

Sources  of  experimental  uncertainties  in  these  results  include  absolute 
energy  calibration  (+8*),  cooling  curve  extrapolations  (+8X),  beam  spot 
area  error  for  conversion  to  fluence  (+5%),  and  variation  in  laser  pulse 
time  profile  for  conversion  to  peak  intensity  (+10$).  Total  estimated 
relative  uncertainties  are  +9%,  with  absolute  uncertainties  of  +15$. 

Figure  35  shows  a comparison  of  the  thermal  coupling  coefficient  for 
beam  spcts  differing  in  diameter  by  a factor  of  two.  The  circles  indicate 
dafca  obtained  with  beams  of  0.46  cm  diameter  and  are  repeated  from  figure 
31.  The  squares  give  the  results  of  similar  coupling  measurements  on  2.22 
cm  diameter  targets  with  beam  spots  of  0.92  c.tu  diameter.  The  data  indicate 
a slightly  lower  value  of  threshold  fluence  as  well  as  smaller  peak  values 
of  thermal  coupling  for  the  larger  beam  spot  size,  but  more  measurements 
are  needed  to  show  this  in  detail. 

c.  Effect  of  Ambient  Pressure 

The  dependence  of  the  thermal  coupling  coefficient  a on  ambient 
pressure  at  the  target  is  shown  in  figure  36  fo*'  aluminum.  The  measurements 
with  beam  spot  diameter  0.46  cm  extend  from  1 to  10  atm.  Tna  data  are 
all  for  laser  pulse  energies  near  130  J except  the  76  Torr  point  which  was 
for  3 112  J.  The  general  trend  of  decreasing  a with  increasing  pressure 
is  associated  with  the  increasing  density  and  shielding  effect  of  the  sur- 
face plasma  formed  at  these  flue.nces.  Similar  behavior  is  seen  for  incite 
mass  removal  over  this  same  range  of  ambient  pressure  (see  Section  V). 
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TABLE  6 

THERMAL  COUPLING  FLUENCE  THRESHOLDS  COMPARED  TO  f1  =0.5  PINHOLE  THRESHOLDS 
FROM  TABLE  4. 


Large-spot 

f\  = 0.5 

pinhole 

threshold 

incident 

fluence 

from 

Table  4 
(0/cm2) 

Thermal 

coupling 

half-step- 

height 

incident 

fluence 

Constant 
energy 
deposi ted 
in  the 
a E=const 
region 

(J/cm2) 

(0) 

Aluminum 

120 

120 

6.1 

Titanium 

84 

95 

10.4 

Stainless 

96 

95 

7.6 

Nickel 

118 

120 

4.8 
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Figure  36.  Dependence  of  Aluminum  Thermal  Coupling  on 
Ambient  Pressure  at  the  Target 


d.  Effect  of  Pulse  Length 

The  effect  of  HF  laser  pulse  length  on  total  thermal  coupling  coeffi- 
cients was  studied  for  aluminum  and  titanium  targets  in  1 atm  pressure  air 
for  the  0.46  cm  diameter  focal  spot.  The  two  temporal  profiles  used  are 
designated  "short  pulse,"  referring  to  the  standard  HF  pulse  with  4 usee 
FWHM  and  about  7 usee  total  length,  and  "long  pulse,"  with  7 ysec  FWHM 
and  11  usee  total  length.  Aluminum  and  titanium  results  are  given  in 
figures  37  and  38.  The  short-pulse  results,  plotted  as  circles,  are  the 
subset  of  those  plotted  In  figures  31  and  34  which  correspond  to  a constant 
pulse  length.  The  long-pulse  data  are  given  as  triangles. 

The  data  above  35  J,  and  hence  well  above  plasma  threshold,  have  been 
least-squares  fitted  to  the  power  law 

a = a Efb  (31) 

The  resulting  curves  are  shown  in  the  plots.  The  short-pulse  thermal 
coupling  gives  b - 0.93  and  1.04  for  aluminum  and  titanium,  corroborating 
the  b = 1 fits  used  In  figures  31  and  34.  There  is  no  apparent  threshold 
shift  due  to  changing  pulse  length.  For  both  materials  the  long-pulse 
coupling  above  threshold  is  larger  in  magnitude  than  short-pulse  coupling, 
and  also  shows  less  shielding  with  Increasing  energy.  The  exponents  of 
Incident  energy  b * 0.84  for  aluminum  and  0.7  for  titanium  indicate  that, 
unlike  short-pulse  behavior,  the  long-pulse  absolute  energy  deposited  in 
the  target  continues  to  Increase  with  incident  energy  even  above  the  plasma 
threshold.  The  total  amount  of  data  taken  here  does  not  justify  the  con- 
clusion that  longer  pulses  are  preferred,  however.  This  Is  true  especially 
because  we  have  only  measured  the  total  ?s  opposed  to  the  local  coupling. 

2.  SPATIAL  DISTRIBUTION  OF  THERMAL  COUPLING 

A topic  of  current  Interest  Is  the  spatial  distribution  of  thermal  energy 
deposited  In  a metal  target  by  laser  beam  interaction.^  Beginning  with  a 
deposition  profile  which  corresponds  to  the  Incident  beam  distribution  at 
low  Intensities,  increasing  energies  lead  eventually  to  plasma  Ignition  and 
the  surface- heating  mechanisms  associated  with  a radiating  high-temperature 
gas  which  Is  spreading  radially  as  well  as  propagating  axially  up  the 
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Incident  beam.  Thermal  deposition  measurements  at  beam  energies  above  LSA 
wave  thresholds  have  shown  heating  profiles  with  widths  greater  than  the 
Incident  beam  width  as  discussed  below  and  in  Section  VI. 

a.  Target  Size  Effects 

Total  thermal  coupling  coefficients  were  measured  as  described  in 
Section  IV  for  targets  of  three  different  diameters  in  order  to  obtain 
information  on  the  radial  extent  of  the  target  heating  profile.  Aluminum 
targets  with  diameters  d]  * 0.79  cm,  d2  = 1.59  cm,  and  d*  ~ 2.22  cm  were 
exposed  to  HF  laser  oeam  pulses  near  80  J with  a oeam  focal  spot  diameter 
of  0.46  cm.  As  seen  in  figure  31,  these  conditions  are  well  above  plasma 
threshold,  lying  in  the  aE  * constant  region.  For  constant  incident 
energy  E^,  coupling  coefficients  a were  measured  to  form  the  ratios  ai/ct2  * 
Wj/V/2  and  (*2/03  = W2/W3  where  Wj  is  the  energy  deposited  in  target  j,  and 
the  subscripts  are  the  target-size  Indices. 

For  purposes  of  defining  a time-averaged  effective  heating  width 
parameter,  an  axially  symmetric  gaussian  with  1/e  radius  6 was  assumed 
for  the  thermal  deposition  profile.  (Also  see  Section  VI.)  An  effective 
gaussian  heating  radius  6 can  be  obtained  from  any  two  target  sizes  by 
solving  for  <5  by  the  relation 

1 - exp  { -r2  / 52  ) W1  (32) 

1 - exp  ( -r2  / 62  ) W2 

where  the  radii  are  known  and  the  ratio  Wj/W2  is  measured. 

The  results  of  table  7 show  that  the  time-averaged  heating  profile 
Is  about  twice  the  diameter  of  the  beam  spot.  The  half -maximum  radius 
r (heating)  * In  2 6.  The  variation  of  the  effective  gaussian  radius 
6 with  the  size  of  the  target  pairs  is  Interpreted  as  reflecting  the 
inadequacy  of  the  single-gaussian  profile  assumption. 
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Targets  2 and  3 (large  pair)  give  a gausslan  width  5 « 0.73;  targets 
1 and  2 (small  pair),  6 ■ 0.44  cm.  This  Indicates  an  actual  deposition 
profile  which  has  a high-level  central  component  and  a lower-1 eve^  broader 
component.  This  particular  departure  from  the  single  gaussian  $r/ape  Is 
also  Indicated  by  time-dependent  temperature  measurements  discussed  below. 


b.  Target  Temperature  as  a Function  of  Time  a/id  Radius 
A second  method  of  obtaining  the  spatial  distribution  of  thermal  energy 
deposited  In  a metal  target  Is  by  temperature  measurement.^  The  aluminum 
targets  used  are  thin  so  that  axial  therw/.l  diffusion  times  are  small 
(typically  > 1 msec),  and  large  enough  1 r,  the  radial  direction  to  be 
effectively  Infinite  as  far  as  thermal  behavior  near  the  center  is  concerned, 
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- 0.40  cm. 


EFFECTIVE  GAUSSIAN  HEATING  RADIUS  5 IMPLIED  BY  TARGETS  OF  RADII  r_ 
r2  ■ 0.80,  AND  r3  ■ 1.11  cm.  FOR  INCIDENT  HF  tf  'ER  FLUENCE  e --  480  J/cm2 
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An  established  method  of  obtaining  the  initial  energy  deposition  profile 
T (r,0)  is  to  measure  the  temperature  at  the  target  center  as  a function  of 
time  and  infer  the  initial  distribution  by  a Laplace  transform.  • ' Character- 
istics of  initial  amplitude  and  width  are  then  obtained  for  an  assumed 
standard  profile,  which  is  typically  a gaussian. 


Here  this  method  has  been  generalized  by  Including  time-dependent  measure 
ments  of  temperature  at  points  away  from  the  target  center.  Interpretations 
are  again  carried  out  in  terms  of  the  gaussian  profile  model.  In  principle 
the  time-dependent  temperature  at  any  radial  distance  from  the  target  center 
may  be  used  to  infer  the  Initial  temperature  distribution: 


T(rQ,  t)  -»■  T(r,  0) 


The  time-dependent  temperature  for  a gaussian  initial  profile  is: 


T.t  r.r#2  1 

To  L 


(33) 


Where  * is  diffusivlty,  Tn  * T(0,  0),  and  t specifies  the  gaussian 

1 ° 

- width  6 at  t * 0: 


= 4 


KT 


(34) 


■ Equation  33  is  plotted  in  figure  39  for  four  values  of  rQ. 


The  maximum  ^perature  for  each  of  these  curves  (for  rQ  - 6)  occurs 
at: 


max 


o 

TT 


(35) 


Experimentally  the  value  of  t|nax  Is  readily  obtained  and  is  Insensitive  to: 
(1)  the  thermocouple  time  constant,  which  must  be  taken  Into  account  at 
early  times,  and  (2)  baseline  zero  shifts  to  which  late-time  data  are  most 
sensitive.  Since  rQ  is  determined  by  thermocouple  location  and  * is  a 
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Figure  39.  Normalized  Temperature  as  a Function  of  Radius 
and  Time,  Applicable  for  a Gaussian  Initial 
Temperature  Profile  In  a Thin  Infinite  Sheet  of 
Aluminum 


function  of  the  material,  an  approximate  value  of  t can  be  obtained  from  each 
curve  of  measured  T(rQ  * 0,  t)  simply  from  its  t . This  maximum  is  expected 
to  occur  only  if  rQ2  > 4kt.  A better  value  of  t can  be  obtained  by  fitting 
the  entire  curve  T(r  , t). 

Measured  temperature  T(rQ,  t)  for  aluminum  are  shown  in  figure  40  for 
laser  pulses  near  140  J.  The  general  trends  expected  from  figure  39  are 
observed,  with  a thermocouple  time  constant  near  10  msec  affecting  the 
early  time  data,  particularly  for  rQ  = 0.  The  value  of  x which  gives  the 
initial  profile  width  is  best  obtained  from  equation  35  using  a value  of 
rQ  for  which  t is  small.  (For  large  values  of  t , one  subtracts  two 
large  numbers.)  A value  of  6 = 0.29  cm  gives  x = 40  msec,  and  predicts 
Snax  ^7,  and  256  msec  for  rQ  - 0.4,  0.6,  and  0.8  cm.  This  agrees 

reasonably  well  with  the  values  of  t^  in  figure  40.  (The  focal  spot 
radius  is  0.23  err,.)  The  fact  that  the  rQ  = 0 data  at  late  times  indicates 
a larger  value  of  x and  larger  6 is  due  to  an  actual  initial  profile  which 
has  a high-level  central  component  as  well  as  a lower-level  broader  compon- 
ent. Therefore,  the  single  gaussian  approximation  is  not  always  adequate. 

c.  Central  Temperature  as  a Function  of  Time 

Setting  rQ  = 0 in  equation  33  gives  the  time  dependence  of  temperature 
at  the  center  of  the  thermal  deposition  profile.  This  "central"  temperature 
was  observed  as  a function  of  time  for  200  msec  for  a 0.46  cm  diameter  beam 
spot  on  aluminum  samples.  Table  8 shows  the  result  of  a two-parameter  fit 
to  <«> 

T-kl 

' ” t+  t (36) 

over  the  time  range  20  to  200  msec  for  three  incident  energies:  8,  54,  and 
145  J.  The  parameters  x and  6 are  related  by  equation  34.  TQ  is  the 
initial  temperature  rise,  Ed  is  the  total  deposited  energy,  and  ac  is  the 
corresponding  thermal  coupling  coefficient.  For  comparison,  total  measured 
thermal  coupling  coefficients  as  obtained  from  figure  31  are  tabulated  as 
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TABLE  8 


HESULTS  OF  TWO-PARAMETER  (T  , t) 
MEASURED  AT  r = 0.  u 


FITS  OF  THE  TIME -DEPENDENT  TEMPERATURE 


Ei 

(J) 

T 

(msec) 

5 

(cm) 

To 

(°K) 

Ea 

(0) 

ac 

at 

8.1 

180 

0.61 

4.6 

0.32 

0.039 

0.051 

54 

170 

0.61 

118 

7.7 

0.14 

0.113 

145 

270 

0.76 

89 

9.3 

0.064 

0.042 
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SECTION  V 

MASS  LOSS  AND  TARGET  IMPULSE  DUE  TO  HF  AND  DF  RADIATION 
1 . MASS  LOSS 

Mass  removal  of  two  ablators  by  pulsed  HF  radiation  was  investigated 
for  large  spot  conditions  at  1 atm  air  pressure.  Lucite  and  Teflon  samples 
were  exposed  at  the  focal  plane  of  a 76  cm  focal  length  mirror.  The  focal 
spot  was  described  in  Section  II  and  had  an  area  inside  the  half-maximum 

p 

intensity  contour  of  0.69  cm  . 

The  energy  range  available  allowed  the  study  of  mass  removal  for  total 
pulse  fluences  from  0 to  250  J/cm^  with  peak  intensities  up  to  75  MW/cm*' . 

In  this  region,  inspection  of  exposed  Lucite  and  Teflon  targets  showed 
distinctly  different  types  of  target  response.  Lucite  mass  removal  was 
observed  as  a shallow  crater  with  approximately  the  shape  of  figure  8 with 
effects  confined  to  the  region  near  the  surface.  Teflon,  on  the  other  hand, 
transmits  most  of  the  HF  radiation  to  a significant  depth.  Damage  in  this 
case  appeared  as  a blackening  of  the  white  Teflon  samples  throughout  the 
0.16  cm  thickness.  The  only  apparent  mechanism  for  mass  loss  is  localized 
small-scale  blistering  on  the  target  front  surface  and,  to  a smaller  extent, 
on  the  rear  surface. 

Mass  loss  curves  are  shown  in  figure  41  as  a function  of  incident 
energy  for  both  materials.  The  mass  loss  of  Teflon  is  zero  on  this  scale. 
The  solid  circles  show  the  specific  energy  requi;  ed  to  remove  Lucite,  in 
kvl/g. 


The  dependence  of  Lucite  mass  loss  on  the  ambient  pressure  p was 

investigated  in  the  vacuum  vessel  of  figure  15.  The  mass  loss  was  measured 

-4 

at  four  pressures  from  1 to  10  atm  corresponding  to  altitudes  of  about 
0,  16,  31,  and  49  km.  The  open  circles  of  figure  42  are  for  85-J  incident 
laser  pulses.  (The  point  at  7.6  Torr  Is  an  interpolation  between  91  J and 
59  J points,  which  are  indicated  by  rectangular  symbols.  Thr  measured 
energies  for  the  other  points  are  85+1  J.)  The  decrease  in  the  density 
and  the  opacity  of  laser-formed  surface  plasma  results  in  increased  coupling 
and  mass  loss  at  reduced  pressures. 
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2.  IMPULSE 

The  0.94-cm  diameter  HF  laser  focal  spot  was  also  used  for  measurements 
of  the  momentum  imparted  to  several  materials  in  air.  A linear  velocity 
transducer  (LVT)  provided  a continuous  record  of  the  instantaneous  velocity 
of  the  target  and  supporting  plunger  through  about  0.7  cm  of  travel  occurring 
over  a few  milliseconds  after  the  laser  shot.  The  voltage  output  of  the  LVT 
was  recorded  as  an  oscillograph.  Because  the  distance  of  travel  was  well 
defined  and  readily  measured,  the  LVT  calibration  could  be  monitored  for 
each  shot. 

Initial  impulse  measurements  were  made  to  determine  the  effect  of  target 
geometry.  All  targets  were  2.54  cm  in  diameter.  Aluminum  targets  were 
exposed  to  pulses  between  150  and  180  J for  two  target  geometries:  (1)  The 
target  was  mounted  on  the  LVT  moving  core  by  a cylindrical  post  to  the 
center  of  the  target  rear  surface  so  that  the  target  edge  was  exposed. 

In  this  case  pressure  waves  occurring  on  the  front  surface  flowed  around 
the  target  edge,  and  decreased  the  net  impulse  delivered  to  the  target 
(2)  The  target  was  mounted  as  in  (1)  except  for  the  addition  of  an  edge 
shield  consisting  of  a large  planar  surface  with  a 2.54  cm  hole  through 
which  the  target  could  move.  In  this  second  geometry,  the  target  was 
initially  flush  with  the  shield  surface,  forming  an  almost  uninterrupted 
plane.  Radial  clearance  was  typically  0.015  cm.  The  shield  plane  was  0.1 
cm  thick.  The  impulse  measured  with  the  shield  in  place  was  typically 
larger  by  a factor  of  1.4  than  that  observed  without  the  shield.  The 
measurements  reported  here  were  all  taken  with  the  shield  in  place. 

The  specific  impulse  delivered  to  aluminum  6061  is  shown  in  figure  43 
as  a function  of  incident  energy  for  both  HF  and  DF  radiation.  The  DF  and  HF 
thresholds  are  both  near  50  J,  but  the  DF  curve  shows  a step  increase  while 
the  HF-produced  impulse  increases  gradually  up  to  about  120  J.  These  Impulse 
curves  are  similar  to  the  OF  and  HF  threshold  behavior  shown  in  figure  18 
which  showed  the  beam  plasma  cut-off  for  the  same  target  material  and  the 
same  beam  spot  size  and  configuration.  For  DF,  both  plasma  cut-off  and 
Impulse  show  a steep  threshold  near  50  J (70  J/cm^).  For  HF.  both  plasma 
cut-off  and  Impulse  show  an  initial  threshold  at  about  the  same  energy,  but 
a gradual  variation  with  energy  over  the  range  50  to  120  J. 
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Figure  43. 


Specific  Impulse  Imparted  to  Aluminum  by 
KF  and  DF  Pulsed  Laser  Radiation 


Less  detailed  momentum  measurements  were  made  for  HF  radiation  incident 
on  stainless  steel,  titanium,  and  Lucite.  These  data  are  plotted  in  figure 
44  with  the  aluminum  HF  data  of  figure  43  as  a background  for  comparison. 

The  four  materials  show  similar  behavior  in  the  energy  range  studied  except 
for  Lucitc  where  vaporization  contributes  to  the  impulse  delivery.  Except 
for  titanium,  the  approximate  impulse  thresholds  observed  correlate  well 
with  the  plasma  threshold  observed  in  figures  21,  22,  and  26. 

3.  IMPULSE  FROM  CONFINED  PLASMA 

Values  of  I/E  about  one  order  of  magnitude  greater  than  the  largest 
values  shown  in  figure  44  have  been  observed  for  a target  geometry  which 
results  in  confinement  of  the  target-produced  plasma.  The  resulting  large 
values  of  delivered  impulse  can  be  expected  wherever  an  absorbing  surface 
is  covered  by  an  IR-transparent  material  having  significant  mass. 

A 0.16  cm  thick  teflon  wafer  was  attached  to  an  aluminum  disc  (both 

2.54  cm  diam)  with  an  opaque  adhesive.  An  HF  laser  pulse  of  143  J delivered 

to  the  teflon  side  of  the  assembly  imparted  a specific  impulse  I/E  of  75 

dyne-sec/J  to  the  aluminum  backing  disc.  This  event  shattered  the  teflon 

disc  and  deformed  both  the  al.iminum  backing  disc  and  the  LVT  plunger.  The 

velocity  of  the  backward-moving  assembly  reached  1.9  x 10  cm/sec.  If 

uniform  acceleration  for  30  psec  is  assumed  an  acceleration  of  6 x 10^  g's 

is  obtained.  Distributing  the  force  over  the  entire  2.54  cm  disc  diameter 

gives  an  average  pressure  during  that  30  usee  of  70  atm;  if  the  force  is 

2 

applied  only  over  the  beam  spot  area  of  0.69  cm  the  average  pressure  is 
730  atm. 
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SECTION  VI 

LASER  INTERACTION  AT  10.6  MICRONS  WAVELFNGTH 


Although  the  primary  emphasis  In  the  experiments  of  this  report  was  or* 
Interactions  at  chemical  laser  wavelengths,  a number  of  new  experiments  were 
carried  out  at  the  longer  COg  laser  wavelength  in  order  to  determine  some 
aspects  of  wavelength  scaling  and  on  mechanisms  of  enhanced  coupling.  From 
experiments  reported  here,  using  different  size  targets,  we  have  determined 

v . 

the  dependence  of  the  plasma  spreading  on  delivered  laser  energy.  These 
results  were  found  to  be  In  approximate  agreement  with  new  experiments  using 
a transverse  air  flow  facility.  With  a high  speed  wind,  the  plasma  is 
blown  partly  away  from  the  desired  targsi  area  and  the  net  coupling  is 
reduced.  / - 

We  also  attempted  experiments  using  moTten  targets  to  determine  those 
conditions  under  which  the  pressure  from  LSD  wav*^  would  blow  out  the  molten 
material. 


1.  EXPERIMENTAL  APPARATUS 


These  measurements  were  made  using  the  Boeing  Double  Pulse  CO^  lasers 
whose  focal  spot  Intensity  profiles  have  been  measured.  The  focal 
spot  has  a 0.16  cm  radius  for  a 37.5  cm  focal  length.  Roth  lasers  have  an 
approximate  6 nsec  full -width  half  maximum  temporal  width,  with  peak  power 
at  1.5  psec  and  a long  tall.  Laser  energy  in  each  beam  can  reach 

15  J with  an  arbitrary  delay  between  the  two  pulses.  In  the  work  reported 
here,  however,  only  single  pulses  or  simultaneous  double  pulses  were  used 
with  each  beam  Incident  at  9°  to  their  common  axis  which  was,  in  turn, 
perpendicular  to  the  surface  of  the  target.  The  operation  of  the  lasers 
was  always  unchanged  as  far  as  discharge  voltage  or  gas  mixture  is  concerned 
and  energies  were  varied  by  the  Insertion  of  attenuators.  Laser  power  and 
energy  were  monitored  by  beam  sampling  diagnostics  which  viewed  the  beams 
after  attenuation.  In  the  focal  plane  there  were  no  hot  spots. 


15.  Hall,  R.B.,  Maher  W.E.,  and  Wei.  P.S.P.,  "An  Investigation  of  Laser- 
Supported  Detonation  Waves."  AFWL-TR-73-28,  June  1973. 

18.  Hall,  R.B.,  Maher,  M.E.,  and  Nelson,  D.J.,  "Double-Pulse  Laser  Inter- 
action Experleents,"  AFMl-TR-73-296,  July  1974. 
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A transverse  air  flow  facility  was  constructed  during  this  study.  Its 
flow  field  was  examined  so  that  laser  target  interaction  experiments  could  be 
performed  in  a well  defined  uniform  transverse  air  flow.  A diagram  of 
the  subsonic  transverse  flow  apparatus  is  shown  in  figure  45.  The  air  is 
stored  in  a 1000  liter  vessel  at  6 atmospheres.  A 10:1  area  ratio  is  used 
for  the  long  profile  nozzle  which  leads  to  the  open  test  volume.  Beyond  that 
a catcher  and  diffuser  section  recovers  the  air  and  sends  it  through  a filter. 
The  high  flow  regulator  valve  can  be  adjusted  for  plenum  pressure  and  hence 
average  velocity,  leading  to  a limited  time  available  for  nearly  uniform  flow. 
After  the  on-off  valve  is  opened,  approximately  five  seconds  are  required  for 
the  flow  to  stabilize.  Over  the  next  ten  seconds  the  average  velocity  is 
nearly  constant,  varying  less  than  T'-.  The  exit  Gf  the  nozzle  is  2 x 2 cm 
and  the  open  region  length  is  4 cm.  Figure  46  is  a ohotograph  of  the  entire 
facility  with  the  storage  tank  in  the  lower  foreground  and  with  the  high 
flow  regulator  valve  in  the  last  leg  in  the  piping  (seen  angling  to  the  upper 
left)  before  the  plenum.  The  experimental  area  is  between  the  nozzle  and 
the  catcher  in  the  upper  right  portion  of  the  picture,  just  to  the  left  of 
the  mirror.  This  a^ea  is  shown  in  an  enlarged  view  in  figure  47.  In  the 
foreground  is  the  large  focusing  mirror,  which  focusses  both  beams  on  a 
position  shown  by  the  bright  spot  in  the  center  of  the  wall  forming  one 
side  of  the  test  volume.  It  is  fastened  to  the  nozzle  outlet  which  is  in 
the  tapering  metal  tube  near  the  right  side  of  the  photograph.  The  flow 
moves  from  left  to  right  through  the  test  volume  as  it  is  shown  in  this 
photograph.  Figure  48  is  a line  drawing  of  the  same  scene  as  figure  47 
with  the  laser  beam  geometry  also  indicated.  These  two  parallel  beams 
enter  the  scene  from  the  right  with  beam  #1  on  top.  After  skirting  the 
test  volume,  nozzle,  and  catcher,  they  reach  the  mirror  which  then  sends 
them  to  the  common  focal  point. 

The  test  volume  was  probed  with  a 0.1  cm  diameter  tube  pointed  upstream 
and  connected  to  a pressure  meter.  Each  datum  point  was  recorded  10  seconds 
after  flow  Initiation  and  determined  the  local  velocity  to  within  one  percent. 
In  addition,  flow  fluctuations  contributed  2 to  3X  uncertainty.  These  data 
were  evaluated  to  find  the  uniformity  conditions  In  the  test  volume  which  are 
stamrlzed  In  figure  49.  Within  the  shaded  areas  In  the  planes  0,  1,  2,  and 
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Figure  45.  Subsonic  Transverse  Flow  Apparatus 


93 


Figure  48.  Lass.-  Beams  and  the  Subsonic  Transverse  Flow 
Experimental  Arrangement 
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WITHIN  oHAOED  AREA  160<u<172  m/ttC,  i.t.,  166  m/*ec +4%. 
NOTE:  T~270°K 


Figure  49  Velocity  Field  of  Subsonic  Flow 
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3 cm  downstream  from  the  nozzle  exit  the  velocity  field  u Is  166  m/sec 
± 4*  at  T * 270°K.  These  measurements  were  made  to  within  0.1  cm  of  the 
target  plane.  The  magnitude  of  the  boundary  layer  <5^  was  calculated  from 


where  u 3 1.82  x 10'4  poise  and  p * 1.22  x 10"3  gm/cm3. 

In  equation  37,  u Is  the  measured  flow  velocity  and  x Is  the  distance 
from  the  point  at  the  nozzle  throat  wbsre  the  boundary  layer  begins.  In 
the  experiment  ■ -ported  here,  x * 5 cm,  indicating  <5^  3 0.034  cm. 

2.  THERMAL  COUPLING  IN  THE  ABSENCE  OF  TRANSVERSE  AIR  FLOW 

Thermal  coupling  of  energy  in  laser  beams  to  targets  Is  dependent  upon 
many  parameters  Including  the  itrlnslc  absorptivity  at  the  laser  wavelength, 
the  formation  of  plasmas  above  the  target  surface,  and  upon  the  subsequent 
plasma  spreading. Several  techniques  have  been  used  in  this  study  to 
distinguish  between  some  of  the  processes  and  to  understand  the  way  plasma 
spreading  takes  place. 

In  these  Investigations,  0.012  cm  diameter  chromel-alumel  thermocouples 
were  spot  welded  together  on  the  target  samples  so  that  the  thermal  junction 
was  at  the  rear  surface.  In  all  the  tests  discussed  In  Section  VI,  the  gross 
thermal  coupling  was  determined  by  measuring  the  equilibrium  temperature  of 
the  target. ^ The  equilibrium  time  for  the  1.59  cm  diameter  targets  used 
here  was  found  to  be  approximately  0.5  sec  for  aluminum.  A 10  second  or 
longer  cooling  record  was  used,  however,  to  give  an  accurate  value  for  the 
equivalent  equilibrium  temperature  rise  at  the  time  of  the  laser  beam 
Irradiation. 

In  figure  50,  thermal  coupling  to  nickel  Is  shown  for  two  focal  lengths. 
The  spot  sizes  art  the  areas  corresponding  to  the  full  width  half -maximum 


17.  Plrrl,  A.  N.,  "Enhanced  Thermal  Coupling",  Second  DOD  High  Energy  Laser 
Conference,  Colorado  Springs,  November  1976. 
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intensity.  The  total  coupling,  as  discussed  above,  had  a maximum  value  of 
0.25  for  both  spot  sites.  This  represents  an  enhancement  of  about  a factor 
of  five  over  the  Intrinsic  absorptivity  of  5X  at  10.6  p.  The  threshold  for 
this  behavior  can  be  seen  to  occur  at  25X  lower  peak  fluence  for  the  larger 
laser  spot  area.  However,  the  monotonlcally  decreasing  absorptivity  keeps 
the  larger  spot  coupling  below  that  of  the  smaller  spot.  The  continued  high 
coupling  at  large  Intensities  arises  because  the  plasma  Is  spreading  out 
across  the  surface  and  heating  the  target  outside  the  laser  spot.  Since 
the  larger  laser  beam  spots  on  these  1.59  cm  diameter  targets  lead  to  somewhat 
smaller  values  of  a,  one  speculates  that  a will  tend  to  be  a constant  for 
large  beam  areas. 

Similar  data  for  titanium  are  shown  in  figure  51.  Again,  the  data 

obtained  with  the  larger  laser  spot  have  a lower  thermal  enhancement  threshold. 

In  addition,  lower  values  of  a throughout  the  enhancement  region  are  obtained. 

These  features  indicate  plasma  spreading  where  a part  of  the  heating  profile 

for  the  larger  laser  spot  extends  beyond  the  1.59  cm  diameter  target.  Data 

(41 

presented  in  a previous  report'  ' for  magnesium  showed  the  same  kind  of 
behavior  for  laser  spot  radii  of  0.16  and  0.35  cm  on  the  same  size  targets 
(1.59  cm  diameter). 

In  order  to  determine  the  magnitude  of  the  plasma  spreading,  the  thermal 
coupling  was  measured  for  target  disks  of  different  dlaiwaters.  A wide  heating 
profile  will  couple  a larger  thermal  energy  into  a larger  diameter  target 
assuming  all  other  conditions  are  constant.  The  approach  used  here  was  to 
obtain  data  for  such  a variation  of  target  diameter  and  then  to  determine  a 
single  parameter  for  the  heating  profile. 

For  these  measurements,  targets  of  aluminum  and  titanium  with  diameters 
of  0.80  and  1.59  cm  were  used.  All  of  the  measurements  were  otherwise  similar. 
The  focal  spot  of  0.16  cm  radius  was  produced  by  a 37.5  cm  focal  length 
mirror.  Each  target  had  a single  centered  chrome! -alumel  thermocouple  and 
was  used  as  received  except  for  a methanol  wash.  The  thermal  enhancement 
threshold  was  about  4 J for  aluminum  and  about  6 J for  titanium  for  this  spot 
size. 
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Since  the  laser  energy  varies  somewhat  from  shot  to  shot,  it  was  of 
interest  to  fit  analytic  expressions  “>  v.he  gross  thermal  coupling  variation 
with  laser  energy.  Thus,  data  taken  w.ch  different  size  targets  could  be 
easily  compared.  For  1.59  cm  diameter  targets: 

o » .29  exp  (-.078  E^  (38) 

(for  aluminum) 


cx  - 1.83  E^*93  (39) 

(for  titanium) 

In  equations  38  and  39,  E^  Is  the  total  Incident  laser  energy  while  the  peak 

fluance  e Is 

max 

emax  a ^*043  joules/cm2  (40) 


Figure  52  gives  plots  of  the  ratio  of  the  thermal  coupling  coefficient 
obtained  with  the  two  sizes  of  targets.  This  ratio  a-j  5g(E.|)/a0  79^1 )»  is 
unity  below  the  threshold  for  plasma  Interactions.  The  titanium  data  are 
represented  by  x symbols,  and  the  aluminum  data  by  circles.  Both  data  sets 
show  that  plasma  spreading  Increases  with  Incident  energy  In  an  approximately 
linear  manner,  after  a step  appears  at  their  respective  thresholds.  A value 
cf  this  ratio  greater  than  unity  means  that  thermal  energy  Is  being  deposited 
In  the  area  outside  of  the  smaller  target.  In  any  case,  since  the  focal 
spot  area  is  much  smaller  than  even  the  smallest  target,  It  Is  clear  that 
a considerable  plasma  spreading  occurs.  In  addition,  since  the  ratio 
Increases  with  Increasing  energy,  the  spreading  also  Increases  with  Increasing 
Incident  energy. 
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Figure  52.  Edge  Effect  In  Enhanced  Thermal  Coupling  for  Two  Target  Materials  i 

i 
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It  Is  assumed  that  the  total  abscrbed  energy  fluence  e(~)  has  the 
Gaussian  form 


e(r)  -A  exp  (-r2/fi2)  joules/cm2 


(41) 


where  A and  6 are  constants  for  a given  laser  energy.  Using  the  data  of 
figure  52,  one  can  then  determine  6 as  a function  of  the  incident  laser 
energy. 

In  *1gure  53,  6 Is  shown  as  a function  of  Incident  energy  above  the 
threshold  for  titanium  and  Indicates  a nearly  linear  dependence  of  6 on 
energy.  Even  at  the  threshold  energy,  6 has  a veiue  three  times  the  laser 
beam  spot  radius  which  is  coincidentally  nearly  the  same  as  the  small  target 
radius.  When  the  Incident  energy  triples,  6 doubles  and  the  coupling  area 
Increases  by  a factor  of  four. 

In  figure  54,  aluminum  Is  the  target.  The  results  are  almost  Identical 
to  those  found  for  titanium.  The  scatter  in  the  data  is  larger  for  aluminum 
than  It  was  for  titanium,  but  the  same  trend  Is  unmistakable  as  the  coupling 
area  Increases  by  a factor  of  four  over  the  range  of  Incident  energies 
Investigated  (4  to  20  J). 

In  an  earlier  experiment,  a Laplace  inversion  of  the  central  time 
dependent  temperatures  was  used  to  determine  at  absorbed  energy  density 
profile  width. ^ These  results  showed  a best  estimate  for  6 of  0.44  cm 
when  9.3  j was  Incident  on  an  aluminum  target  with  a focal  spot  radius  at 
half  maximum  Intensify  of  0.078  cm.  Ignoring,  for  the  moment,  the  difference 
In  focal  spot  size,  this  0.44  cm  value  for  6 Is  nearly  the  same  as  the 
results  shown  In  figures  51  and  52  where  the  average  6 for  the  two  target 
diameters  Is  0.51  cm  for  9.3  0. 
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Figure  54.  Heating  Profile  Gaussian  Radius  6 From 
Aluminum  Target  Diameter  Variation 


3.  THERMAL  COUPLING  IN  THE  PRESENCE  OF  AIR  FLOW 

In  the  experiments  reported  in  this  section,  the  transverse-flow  wind 
tunnel  was  used  to  simulate  a fast-moving  target.  In  addition,  it  was  hoped 
that  the  experimental  results  would  aid  in  understanding  the  mechanisms  of 
enhanced  coupling. 

In  these  experiments  using  the  transverse  air  flow,  however,  adjustments 
had  to  be  made  in  the  experimental  procedures  compared  with  the  no  flow  case. 
There  is  strong  convective  cooling  of  the  target  by  the  air  flow  which  had 
a temperature  of  270°K.  This  flow  was  found  to  cool  the  adjacent  targets 
at  rates  of  -0.6°K  sec"^  by  the  time  5 seconds  has  elapsed.  Most  of  the  data 
taken  below  utilized  an  electrical  compensation  circuit  to  avoid  a reduction 
of  temperature  measuring  sensitivity.  With  the  chart  recorder  dynamic 
voltage  range  available  and  the  sensitivity  setting  normally  used  in  no  flow 
thermal  coupling,  only  H seconds  of  flow  would  cover  the  total  available 
voltage  range.  Of  this,  the  first  5 seconds  during  flow  stabilization  wasn't 
useable.  A ramp  voltage  derived  from  an  oscilloscope  triggered  when  the  flow 
commenced  was  placed  in  series  with  the  input  to  the  chart  recorder  so  that 
after  a predetermined  time,  usually  5 seconds,  the  signal  to  the  recorder 
was  not  changing  rapidly.  Then  a laser- interact! on  thermal-coupling  temperature 
excursion  could  be  easily  triggered  and  analyzed. 

Again,  a single  target-cantered  chromel-alu  lei  thermocouple  was  used.  A 
Plexiglas  holder  gripped  each  target  by  the  edge  and  positioned  it  so  that 
it  was  flush  with  the  holder  surface  which  formed  one  side  of  the  test  volume. 
Every  datum  was  from  a fresh  sample  in  an  as -received  condition  except  for 
a methanol  wash.  Every  test  began  with  the  target  at  room  temperature. 

Maximum  equilibrium  temperature  excursions  due  to  laser  thermal  coupling 
were  under  10  K.  Where  two  simultaneous  laser  pulses  were  used  to  Increase 
the  laser  energy,  an  oscilloscope  monitored  detectors  for  both  beams  and  a 
photographic  record  of  the  oscilloscope  traces  was  used  to  assure  that  both 
pulses  were  simultaneous  to  within  1 usee . 


Two  target  sizes  were  used,  and  flew  could  be  present  or  absent.  In 

figure  55,  data  for  a titanium  target  of  diameter  0.81  cm  are  shown  for  a 

flow  velocity  of  166  msec"^  and  for  the  case  of  no  flow.  It  Is  apparent 

that  the  presence  of  flow  makes  no  appreciable  difference  In  the  thermal 

coupling  measured  at  any  incident  energy  in  the  range  6 to  20  J.  The  thermal 

-2  5 

coupling  decreases  as  E^  * and  leads  to  values  of  a near  0.03  at  the 
highest  energies.  This  0.03  value  is  far  below  the  intrinsic  10.6  p 
absorptivity  of  0.14  for  titanium.  Even  at  the  plasma  threshold,  the  maximum 
coupling  Is  below  the  0.34  peak  values  seen  in  figure  51. 

Using  1.59  cm  diameter  targets  and  for  the  same  laser  beam  focal  spot 
size,  a large  change  In  thermal  coupling  occurred  when  the  transverse  flow 
was  on.  These  data  are  shown  in  figure  56.  The  no-flow  data  are  the  same 
as  those  shown  in  figure  49  and  fall  off  nearly  as  E-”"*.  When  the  transverse 
flow  Is  added,  the  thermal  coupling  values  fall  by  60%.  This  is  In  great 
contrast  to  the  0.81  cm  diameter  target  data  and  occurs  because  of  the 
relatively  large  values  of  6 (the  effective  plasma  coupling  radius).  When 
the  plasma  radius  is  either  very  small  or  large  compared  to  the  target  size, 
the  external  flow  moves  the  plasma  with  respect  to  the  target  during  the 
coupling  time  tc  and  the  same  amount  of  coupling  occurs.  Strictly  speaking, 
this  occurs  when 


ut  « greater  of  ^R,  6) 


(42) 


The  plasma  coupling  time  Is  approximately  equal  to  the  laser  pulse  time. 

For  (R/6)  M,  the  external  flow  displaces  the  coupling  plasma  with 
respect  to  the  target  and  a reduction  In  a Is  found. 

A similar  behavior  Is  found  In  tests  with  aluminum  targets  exposed  to 
both  laser  beam  Irradiation  and  transverse  flow.  The  data  for  the  0.81  cm 
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diameter  targets  are  displayed  in  figure  57.  There  is  an  enhancement  region 
above  the  threshold  at  3.5  to  4 J incident  energy  which  yields  values  of  a 
of  about  0.1C.  Such  a thermal  coupling  maximum  value  lies  below  the  maximum 
a value  of  0.33  presented  in  an  earlier  report.  ' In  this  case,  o decreases 
approximately  as  E~^  and  air  flow  causes  a negligible  change  in  coupling. 

For  1.59  cm  diameter  aluminum  targets,  the  variation  In  thermal  coupling  is 
showri  in  figure  58.  Again,  there  is  a greater  reduction  in  a than  found 
with  small  targets.  This  reduction  of  28%  is  only  about  1/2  of  that  seen  for 
titanium.  The  maximum  a values  are  only  slightly  above  0.2,  and  the  nearly 
E“^  behavior  seen  at  the  higher  Incident  energies  does  not  carry  all  the  way 
down  to  threshold. 

The  reduction  In  thermal  coupling  caused  by  transverse  air  flows  occurs  only 
for  the  larger  targets.  This  observation  can  be  explained  using  a simple 
thermal  coupling  model.  This  coupling  model  assumes:  (1)  that  the  thermal 
coupling  radial  profile  is  a constant  "top-hat"  distribution  of  radius  R}J » 

(2)  that  the  external  transverse  air  flow  moves  this  profile  with  respect  to 
the  target  at  the  velocity  of  the  air  flow,  and  (3)  the  thermal  coupling  to 
the  0.81  cm  diameter  targets  is  unchanged  as  found  by  experiment. 

Figure  59  Indicates  the  geometry  of  the  situation.  The  two  shaded 
concentric  areas  represent  the  two  different  target  sizes  with  air  flow 
Incident  from  the  left.  Ry  is  the  radius  of  the  larger  target.  Representative 
plasma  heating  areas  are  designated  by  the  displaced  circles  of  areas  A^, 

Aj,  and  A^.  These  heating  areas  are  all  shown  displaced  such  as  to  be  tangent 
to  the  smaller  target.  Under  these  conditions,  obviously,  no  decrease  in 
power  transfer  to  the  smaller  target  occurs.  Thermal  heating  profiles  of  radii 
R^  below  3 Rt/4  never  cause  a net  decrease  in  power  transfer  to  the  large 
target.  For  larger  heating  radii,  the  Instantaneous  transverse-flow  caused 
coupling  reduction  y Is  shown  in  figure  60.  Since  the  net  energy  reduction 
to  large  targets  Is  an  Integral  of  the  power  transfer,  figure  60  represents 
an  upper  limit  to  the  actual  reduction  that  can  occur. 
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Figure  58.  Trent verte  Air  Flow  Effect  on  Thermal  Coupling 
To  Aluminum  Targets  of  1.59  cm  Diameter 
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POSITION  OF  TOP  HAT  THERMAL  COUPLING  HEAT  INPUT  AREAS 
AFTER  MOVEMENT  DUE  TO  FLOW  OF  DURATION  X/V 


FLOW 


VELOCITY  V 


ALLOWED  MOVEMENT  X 
FOR  R^|  * R.^ 

ALLOWED  MOVEMENT  X 
FOR  RH  - 3/3  Rj 

ALLOWED  MOVEMENT  Xk 
FOR  Rh  - 2 Rt 


Figure  59.  Geometric  Model  for  Floe  Effect  on  Thermal 
Coupling  Btsed  on  Date  from  Two  Target  Sizes 


rh/rt 


Figure  6C.  Flow-Caused  Reduction  In  Energy  Delivered  to  Target  of  Radius  Rj 
By  a Top  Hat  Thenaal  Coupling  Profile  of  Radius  Ru,  Which  Lasts 
a Time  (R«  - 0.5  «t)/V,  Where  Rh  * 0.5  Rt 


When  R^/Hj.  Is  unity,  this  reduction  reaches  a maximum  value  of  31X.  The 
model  results  shown  are  only  Intended  to  show  graphically  the  combined  effects 
of  transverse  flow  and  plasna  spreading  on  measured  thermal  coupling.  It  Is 
evident  that  net  much  reduction  ever  occurs  for  very  small  targets.  For 
maximum  coupling  reduction  on  the  large  target,  the  plasma  must  be  moved  a 
distance  of  0.4  cm.  This  ImDlles  that,  the  coupling  time  t Is  as  long  as 

V 


0 4 

t ■ i 24  microseconds  (43) 

1.66  x 104 


which  Is  long  compared  to  the  time  of  maximum  luminosity  of  observed  laser- 
supported  detonation  waves/1  ^ but  Is  approximately  equal  to  the  length  of 
the  laser  pulse. 


4.  IMAGE  CONVERTER  PHOTOGRAPHY  AND  THRESHOLDS  IN  THE  PRESENCE  OF  AIR  FLOW 
Image  converter  camera  photographs  with  a resolution  of  0.05  cm  were  taken 

of  laser-produced  plasma  events  botn  with  and  without  flow.  These  all  were 
obtained  using  exactly  the  same  setup  ar.d  delay  times,  and  approximately 
same  Incident  energies  for  the  cases  with  and  without  flow.  The  luminous 
phenomena  projected  more  than  1 cm  from  the  surface  when  the  energies  were  in 
the  laser- supported  detonation  wave  regime.  Experimental  apparatus  allowed 
both  absolute  translation  to  0.05  cm  and  angular  changes  to  3°  to  be  observed. 
During  the  approximately  9 usee  when  luminosity  was  sufficient  for  photography, 
no  'hanges  were  ever  seen  during  flow  as  compared  with  no  flow.  During  the 
9 psec  a velocity  of  166  m/sec  would  move  gases  In  the  flow  region  0.15  cm. 

The  incident  Intensity  or  Incident  fluence  required  for  the  Ignition  of 
LSA  waves  on  aluminum  targets  was  Investigated  for  differences  caused  by  the 
presence  of  an  M * 0.5  transverse  air  flow.  A fresh  unexposed  target  was  used 
each  time.  In  as-received  condition  except  for  a methanol  wash.  Two  separate 
series  of  tests  were  completed  and  the  measured  thresholds  were  identical  both 
with  and  without  flow. 

5.  EXPERIMENTS  USING  MELTED  TARGETS 

In  order  to  simulate  a sequence  of  previous  pulses  of  a repetitively-pulsed 

train,  a series  of  experiments  was  carried  out  using  molten  targets  and  single 

pulses  of  CO  “idlatlon.  One  of  the  primary  goals  of  these  experl ements  was  to 

determine  tk  conditions  under  which  the  pressure  from  either  an  LSD  wave  or 

(181 

the  vaporizing  target  Itself  would  blow  away  the  molten  target  material.  v ' 

In  order  to  permit  experiments  with  variations  In  target  temperature;  gallium. 
Indium,  and  tin  targets  were  obtained.  This  gave  us  the  opportunity  to  do  melt- 
ing point  tests  at  303°K,  430°K,  and  505°K.  Mercury  was  also  available  for  room 
temperature  tests.  The  standard  target  was  0.125  cm  thick  and  1.59  cm  in  diameter. 
Each  was  placed  In  a recess  In  a heating  block  attached  to  a thermostatically 
controlled  heater.  A thermocouple  In  the  heating  block  allowed  Its  temperature 
and  temperature  stability  to  be  measured.  The  arrangement  for  the  incident  laser 
beam,  target,  target  block,  and  heater  Is  shown  In  figure  61.  In  this  arrangement 


18.  Fox,  J.  A.,  "A  Method  for  Improving  Continuous  Wave  Laser  Penetration  of 
Metal  Targets,"  Appl.  Phys.  Lett  26,  682,  1975. 
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Figure  6».  Vertical  Arrangement  for  Melted  Target  Studies 


m 


thr  position  of  the  flow  apparatus  is  undisturbed  with  the  melted  target  in 
the  horizontal  plane  which  forms  the  bottom  side  of  the  open  flow  channel.  A 
vertically  Incident  laser  beam  Is  provided  by  a folding  mirror  which  sends  the 
beam  upward  to  the  50.8  cm  focal  length  focussing  mirror.  This  provides  a focal 
spot  0.22  cm  in  radius. 

Without  describing  In  detail  all  the  laser  experiments  on  melted  targets, 

’t  should  be  noted  that  no  gross  removal  of  metal  target  material  was  ever 
(bserved.  The  main  effect,  with  these  small  spot  COg  laser  beams,  was  to  cause 
an  oxidation  of  the  molten  surface.  In  this  oxidized  surface,  a permanent  set 
of  "wrinkles"  remained  as  a result  of  the  formation  of  the  LSD  wave  pressure 
pulse. 

Following  Ignition  of  a laser- supported  detonation  wave  a high  surface 
pressure  occurs  which  has  a strong  radial  gradient.  When  the  target  is  molten, 
this  pressure  pulse  might  be  expected  to  push  material  away  from  the  focal  spot. 
This  was  not  observed  in  the  proceeding  tests  possibly  due  to  surface  tension 
and  viscosity,  so  two  other  series  of  experiments  were  devised  to  look  for  this 
effect.  In  the  first  series,  horizontal  0.81  cm  diameter  targets  were  melted 
In  holders  which  had  holes  completely  through  the  back.  The  high  pressure  pulse 
was  expected  to  push  some  of  the  target  material  out  this  hole.  Without  irradia- 
tion the  surface  tension  prevented  the  melted  target  from  running  out  the  holes  or 
from  doing  anything  other  than  bulging  slightly  downward.  Hole  sizes  were  0.16, 
0.32,  and  0.48  cm  In  diameter.  In  figure  62  the  results  of  these  tests  are  shown. 
In  every  case,  molten  metal  was  extruded  and  remained  hanging  down  In  a drop  of 
tubular  shape.  The  photographs  show  the  results  after  the  target  and  holder  was 
cooled  back  to  room  temperature.  Evldentally,  target  temperatures  were  Insuffi- 
cient to  allow  anything  but  highly  viscous  motion. 

Additional  tests  were  done  with  water  and  mercury  to  look  at  mass  loss  with 
molten  materials  with  lower  vlfttoslty*  In  tests  with  water,  almost  all  the 
water  In  the  holder  cavity  was  blown  out  by  the  pressure  pulse  from  the  laser 
beam-target  Interaction,  however,  similar  tests  with  mercury  removed  0.06  g of 
material,  approximately  if  of  the  target. 
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SECTION  VII 


MELT  REMOVAL  DYNAMICS  DUE  TO  VAPORIZATION 


1 . INTRODUCTION 

It  is  apparent  that  high  intensity  laser  pulses  lead  to  melting  and  vapori- 
zation. At  relatively  low  Intensities  the  molten  layer  can  flow  off  the  surface 
due  to  gravity  or  be  removed  due  to  aerodynamic  forces. 


(19,  20) 


This  report  summarizes  Initial  concepts  and  calculations  on  the  mechanisms 
by  which  a molten  fluid  layer  can  be  pushed  rapidly  outward  by  a local  pressure 
distribution.  In  general,  the  solutions  of  the  combined  continuity  of  mass, 
continuity  of  momentum,  and  continuity  of  energy  equations  are  complex  for  fluid 
flows.  However,  we  have  developed  relatively  simple  approximate  solutions  for 
special  cases  of  interest  for  high  power  laser  coupling.  A laser-melted  surface 
usually  can  be  described  by  a thin  molten  layer  on  top  of  the  underlying  solid 
target.  In  this  case  (of  a thin  molten  layer),  one  can  describe  the  local 
pressure  distribution  as  having  only  radial  gradients  and  being  uniform  along 
the  axial  direction.  As  a model  of  this  situation  the  molten  layer  is  visualized 
as  a slab  of  Incompressible  fluid  with  constant  thickness.  This  slab  Is  bounded 
by  a plane  rigid  surface  on  one  side  and  with  a parabolic  pressure  variation 
applied  to  the  other  side.  The  source  of  the  pressure  distribution  can  be  con- 
sidered to  be  either  a laser- supported  detonation  wave  (LSD  wave)  or  by  the 
existence  of  the  vaporizing  surface  itself.  In  the  first  case,  the  LSD  wave  is 
produced  In  air  and  generally  spreads  out  over  an  area  greater  than  that  of  the 
layer  focal  spot.  In  the  second  case,  where  vaporization  is  important,  the 
pressure  Is  constrained  to  be  within  the  focal  spot.  In  addition,  this  pressure 
Is  easily  produced  In  vacuum  and  may  be  Important  for  laser  space  applications. 


19.  O'Keefe,  0.  D.  and  Johnson,  R.  L.f  "Laser  Melt- Through  Time  Reduction  Due 
to  Aerodynamic  Melt  Removal,"  AIAA  J.  14,  776,  1976. 

Robin,  J.  E. , and  Nordln,  P.,  “Effects  of  Gravitationally  Induced  Melt 
Removal  on  CW  laser  Melt-Through  of  Opaque  Solids,"  Appl.  Phys.Lett.  27. 
593,  1975. 
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With  high  Intensity  CW  or  pulsed  lasers  one  expects  surface  vaporization 
to  exert  enough  reaction  pressure  on  the  molten  surface  to  blow  away  In  the 
radial  direction  some  of  the  molten  surface.  Thus,  although  the  vaporiza- 
tion of  a unit  mass  of  solid  takes  between  10  and  20  times  as  much  energy  as 
Is  required  to  just  melt  the  solid,  the  vaporization  process  itself  can  aid 
In  removing  the  melt. 

The  analysis  given  In  Sections  7.2  and  7.3  describes  the  manipulation  of 
the  conservation  equations  of  mass,  momentum,  and  energy  to  obtain  approximate 
relations  describing  the  dependence  of  the  melt  recession  velocity  vm  and  the 
melt  thickness  s on  the  laser  Irradlance  and  focal  spot  size.  The  surface 
pressure  variation  Is  analyzed  In  self-consistent  terms  using  the  Clauslus- 
Clapeyron  equation.  At  the  higher  Intensities  (as  described  In  Section  7.3) 
the  mass  loss  *n  vapor  form  becomes  Important  and  leads  to  a significant  decrease 
in  the  effective  overall  mass  removal  efficiency. 

The  results  of  the  analyses  of  Sections  7.2  and  7.3  are  shown  in  figures 
63  through  75.  The  specific  targets  studied  Include  aluminum,  titanium,  and 
molybdenum. 

Section  7.5  discusses  the  efficiency  of  removing  a liquid  melt  layer  by  a 
combined  CW  and  pulsed  laser.  In  this  preliminary  analysis,  the  potential  of 
this  technique  appears  quite  limited.  This  Is  due  to  the  large  melt  removal 
times  (for  large  beam  areas)  compared  to  characteristic  times  required  to  melt 
the  target  Itself. 

Figure  63  shows  the  calculated  steady  state  thicknesses  of  a molten  layer 
of  aluminum  as  a function  of  the  absorbed  Intensity  In  the  center  of  the  laser 
beam  with  the  scale  length  of  the  parabolic  pressure  variation  as  a parameter. 

The  curves  in  this  and  following  figures  are  not  continued  Into  the  range  where 
the  melt  thickness  Is  greater  than  a ter.th  of  the  pressure  radius.  The  steady 
state  melt  thickness  decreases  with  Increasing  absorbed  Intensity.  The  curves 


21.  Robin,  J.  E.  and  Nordln,  P.,  "Improved  CW  Laser  Penetration  of  Solids 
Using  a Superimposed  Pulsed  Laser,"  Appl.  Phys.  Lett.  29,  3,  1976. 
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ABSORBED  INTENSITY  IN  CENTER  OF  BEAM  (W/cm2) 
Figure  63.  Steady- State  Melt  Thickness  - AlufRir.um 


transfer  from  a smaller  thickness  curve  dominated  by  me>t  removal  by  radial 
flow  of  the  liquid  to  a larger  thickness  curve  dominated  by  melt  removal  by 
vaporization.  The  range  of  absorbed  Intensities  over  which  the  transfer  occurs 
depend,  on  tne  pressure  radius  with  the  transfer  occurring  at  lower  intensities 
for  larger  pressure  radii. 

Figure  64  shows  the  calculated  steady  state  molting  velocity  of  &lumini-m. 
The  steady  state  melting  velocity  increases  with  increasing  absorbed  inten  ty. 
The  curves  transfer  rom  a larger  velocity  curve  dominated  by  melt  removal  by 
radial  flow  of  the  liquid  to  a smaller  velocity  curve  dominated  by  melt  removal 
by  vaporization.  Figure  65  shows  the  calculated  specific  melt  volume  of 
aluminum,  the  volume  of  melt  removed  from  the  crater  per  joule  of  absorbed 
energy.  The  curves  again  change  from  a larger  specific  volume  curve  dominated 
by  melt  removal  by  radial  flow  of  the  liquid  to  a smaller  specific  volume 
curve  dominated  by  melt  removal  by  vaporization.  All  the  values  are  consider- 
ably below  the  ideal  one  dimensional  value  with  zero  melt  thickness.  For  a 
certain  pressure  radius,  the  specific  melt  volume  curves  show  maxima  with  respect 
to  the  absorbed  intensity.  The  pressure  radius  is  the  radius  of  the  crater . 

The  radius  of  the  molten  region  is  larger  but  only  the  molten  material  out  to 
the  pressure  radius  is  removed  by  the  radial  pressure  gradient. 

Figure  66  shows  that  the  calculated  pressure  radii  are  approximately 
equal  to  the  laser  beam  radii  for  aluminum  for  the  larger  absorbed  intensi- 
ties but  the  pressure  radii  are  smaller  than  the  laser  beam  radii  for  the 
smaller  absorbed  intensities.  With  laser  beam  radius  as  a parameter,  specific 
melt  volume  curves  do  not  show  maxima  but  decrease  monotonically  with  increas- 
ing absorbed  intensity.  Figure  67  shows  the  calculated  specific  impulse  of 
aluminum,  the  force  on  the  surface  per  watt  of  absorbed  power  or  the  impulse 
on  the  surface  per  joule  of  absorbed  energy.  The  specific  impulse  always 
becomes  nearly  asymptotic  to  a certain  value  for  large  absorbed  intensities. 

For  lower  absorbed  intensities,  when  melt  removal  Is  dominated  by  radial  flow 
of  the  liquid,  the  melting  Is  relatively  efficient.  For  higher  absorbed 
Intensities,  when  melt  removal  Is  dominated  by  vaporization,  the  application 
of  force  to  the  surface  is  relatively  efficient. 
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ABSORBED  INTENSITY  IN  CENTER  OF  BEAM  (W/cm'M 
Figure  66.  Ratio  of  Beam  Radius  to  Pressure  Radius  - Aluminum 


Figures  68  through  71  show  the  steady  state  melt  thickness,  melting 
velocity,  specific  melt  volume  and  specific  impulse  for  titanium.  Figures 
72  through  75  show  the  same  for  molybdenum.  These  materials,  aluminum, 
titanium  and  molybdenum,  were  chosen  to  represent  a low  melting  temperature 
metal,  a medium  melting  temperature  metal  with  very  low  thermal  conductivity 
and  a high  melting  temperature  metal,  respectively.  The  steady  state  melt 
tnickness  of  molybdenum  is  less  than  aluminum  but  titanium  is  still  less  because 
of  Its  low  thermal  conductivity.  The  steady  state  melting  velocity  of  aluminum 
Is  greater  than  titanium  which  is  greater  than  molybdenum  when  melt  removal  by 
vaporization  is  dominant  but  sometimes  molybdenum  melts  faster  than  titanium 
for  smaller  absorbed  Intensities.  The  specific  melt  volumes  ar«  in  an  inverse 
order  of  the  melting  temperatures,  aluminum  greater  than  titanium  which  is 
greater  than  molybdenum.  The  asymptotic  value  of  the  specific  impulse  of 
aluminum  is  larger  than  molybdenum  which  Is  larger  than  titanium  because  of 
its  low  thermal  conductivity.  However  titanium  transfers  from  melt  removal 
dominated  by  radial  flow  of  the  liquid  to  melt  removal  dominated  by  vaporization 
at  smaller  values  of  absorbed  Intensity  than  for  aluminum  and  molybdenum  because 
of  the  low  thermal  conductivity  of  titanium.  For  tha  smaller  range  of  absorbed 
intensities  titanium  has  larger  specific  Impulses  than  both  aluminum  and 
molybdenum. 

2.  ONE  DIMENSIONAL  STEADY  STATE  MELTING 

Consider  a one  dimensional  steaoy  state  melting  situation  of  a semi- 
infinite solid  In  which  the  axial  melting  velocity  is  constant  and  the  liquid 
layer  thickness  Is  kept  constant  by  melt  removal  in  the  radial  direction  by  the 
reaction  pressure  gradient  due  to  evaporation  from  the  liquid  surface.  The 
derivation  is  an  extension  of  the  method  used  by  Masters.^  ^ For  a constant 
melting  velocity,  vm,  the  geometry  Is  shown  In  figure  76  where  the  temperature 
T Is  a function  only  of  x where 

* ‘ z - V (44) 


22.  Masters,  J.I.,  "Problems  of  Intense  Surface  Heating  of  a Slab  Accompanied 
by  a Change  of  Phase"  J.  Appl . Phys.  27,  477,  1956 
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100  CENTIMETERS 


ABSORBEO  INTENSITY  IN  CENTER  OF  BEAM  <W/cmz) 

Figure  73.  Steady-State  Melting  Velocity  - Molybdenum 


ABSORBED  INTENSITY  IN  CENTER  OF  BEAM  (W/cm2) 
Figure  75.  Specific  Impulse  - Molybdenum 


The  positive  directions  for  z and  x are  Into  the  solid.  The  solid-liquid 
Interface  Is  at  x B 0 (or  z * vmt)  and  the  liquid-gas  Interface  Is  at  x = -s, 
corresponding  to  a melt  layer  of  thickness  s. 


In  the  solid,  the  heat  equation,  with  constant  thermal  conductivity  K$, 
specific  heat  C , and  density  p (which  Implies  C * C ) Is 


PC  — 
PS  S 3t 


< 

S 3Z2 


(45) 


In  terms  of  the  x variable  equation  45  becomes 


„ r 3T  „ v 32T 
'vmpsCs  3x  Ks 


(46) 


Equation  46  can  be  Integrated  using  appropriate  boundary  conditions  at  x * 0 

and  ».  In  terms  of  the  solid  Initial  temperature  TQ  and  the  melt  temperature 

T , one  obtains 
m 

T ’ To  * <VV  t - ^ 


(47) 


The  heat  flux  at  x ■ 0 Is 


y 3T 
KS  9Z 


VsCs  'Tm  - V 


z ■ vt 
m 


(48) 


For  solid  aluminum,  the  averaged  thermal  properties  are 


ps  « 2.64  gm/cm 


1.05  joules/gm/K 
2.29  watts/cm  °K 


(49) 


for  the  density,  heat  capacity.  the  thermal  conductivity  respectively. 
Thise  thermal  properties  and  those  for  titanium  and  molybdenum  are  shown,  for 
reference.  In  Table  9.  with  * values  take*  largely  from  reference  23. 
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Boiling  Teaperature  2766.8  3591.  4919.  degrees  Kelvin 

Hmi$  of  Vaporization  at  Tb  10.7768  8.8772  6.1409  Kilo  joules/gram 
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Masters'  solves  the  heat  conduction  equation  in  the  liquid  layer 
using  a very  similar  technique  and  determines  the  surface  temperature  T_ 


Tm  + 


aq_ 


VlCL 


(50) 


where  aq  Is  the  absorbed  laser  Irradiance  and  vm  is  again  the  melt  velocity. 
Using  appropriate  boundary  conditions,  one  arrives  at 


'm 


aq  exp  [ - vms/.cL] 

"A  + »,c,'VYo> 


(51) 


In  equation  51,  nm  is  the  latent  heat  of  melting  and  *L  is  the  dlffusivlty 
of  the  liquid  layer.  It  Is  seen  that  the  melt  velocity  decreases  with  increas- 
ing melt  layer  thickness  s.  For  aluminum  ^nm+Cs(Tm-T0)}  -1088.  joules/gm  and 
one  obtains  the  corresponding  zero- thickness  velocity 


* 3.4xl0"4 


aq  cm/sec 


(52) 


s*0 


Fg>  ..ie  above  analysis  to  be  valid  the  melt  layer  is  assumed  to  be  maln- 

(19) 

tained  constant  by  some  external  force  such  as  gravity,  wind  drag  forces  » 
or  evaporation.  Vie  will  now  examine  the  mechanisms  by  which  the  reaction  force 
of  a vaporizing  surface  can  maintain  a "constant"  melt  layer  thickness.  In 
principle,  a similar  analysis  could  be  carried  out  for  the  case  of  external 
surface  pressures  generated  by  laser-supported  absorption  waves^. 


, oruer  v:  „«cermine  relationships  taking  Into  account  the  surface  pres- 
sure gradients  one  must  solve,  self-consistency,  the  equations  of  conserva- 
tion of  energy,  mass,  and  momentum.  For  analytical  simplification.  It  Is 
assumed  that  the  -.*•*  ice  pressure  is  approximated  by 


23.  Toulouklan,  Y.S.,  Powell,  R.W.,  Ho,  C.Y.,  Klemens,  P.S.,  Nlcolaou, 

M.C. , Kirby,  R.K,  Taylor,  R.E.,  and  Desal,  P.D.,  Thermophysical  Properties 
of  Hatter.  Thermal  Conductivity,  Volume  1.  Thermal  Dlffusivlty,  Volume  10, 
and  Thermal  Expansion,  Volume  12,  IFI/Plenum,  1970-1975. 


Po  0 - 


4-> 


surface 
x ■ -s 


where  Rp  Is  the  scale  length  for  the  pressure  variation.  The  maximum  pressure 
PQ  Is  to  be  determined  by  solving  for  the  vaporization  rates.  The  actual 
scale  length  Rp  Is  found  to  be  considerably  smaller  than  the  laser  spot  size 
due  to  the  strong  dependence  of  the  vapor  pressure  on  the  exact  surface  temp- 
erature. Thus,  It  Is  found  that,  for  r<Rp,  the  maximum  variation  in  surface 
temperature  Is  only  3 to  10  percent. 

Thus,  In  the  conservation  of  energy  equation  we  neglect  ~ and  obtain 

or 

’LCL  ( S * **  5T  > - \ J I 

Equation  54  neglects  effects  of  viscosity  and  the  kinetic  energy  of  the 
fluid  in  comparison  to  the  heat  capacity  of  the  liquid.  The  conservation  of 
mass  equation  (with  similar  assumptions)  Is 

r If  <"V>  * ?!*  * 0 < 


Finally  the  momentum  equations  In  both  the  radial  and  axial  directions  are 


3 V 3V  av„ 
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In  equation  56,  y Is  the  viscosity  of  the  molten  target  material.  Using 
the  fluid  boundary  conditions 
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(57) 


vz  {x*0}  - vr  (r*0}  ■ vr  {x=0}  *0 


vz  (x-s)  - vm 


Using  equation  53  one  finds  from  equation  56  that 


v „ -v 

r m 

2 2 

V * Vm  X /s 

z m 


& 

Ipl 


Finally,  the  total  variation  (within  the  molten  layer)  of  pressure  is  given  by 
p.p  /(,._£  \ ♦ _Jt_  (x+s)  ♦ 4-  (x2-s2)  - I 

0 V'  "Sv.  *T  Rp  5 / 

It  may  be  noted  from  equation  59  that  all  the  terms  describing  the  axial 
pressure  variation  are  small,  being  of  Older  s^/Rp.  In  this  case,  one  may 
assume  that  the  pressure  everywhere  is  the  same  as  the  surface  pressure  (see 
equation  53). 

Knowing  the  axial  fluid  velocity  given  by  equation  58,  one  can  determine 
the  variation  of  temperature  with  the  molten  layer.  Setting  Y * gj-»  one  finds 

31  ♦ ^ (1  - 4 ) Y - 0 

ax  «L  7 

Equation  60  can  be  Integrated  twice  and  with  the  application  of  heat  flux 
boundary  conditions  at  x*0  (the  solid-liquid  boundary)  and  x*-s  (the 
liquid-air  boundary)  one  obtains 
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(61) 


T . T + l«a.  N [A,  ] a) 

ad  <xp  [-  | A]  b) 

v” 


In  equation  61 

vs 

A - JS-  (62) 

kL 


Thus,  Me  find  that  the  melt  velocity  vm  obtained  with  fluid  motion  Is  Identical 
to  equation  51  except  for  the  factor  of  2/3  In  the  exponent.  The  function 
N{A,w>  Is  shown  in  figure  77  and  Is  defined  by 

w 

N {A,w}  - exp  exp  {A(W  - | W2 3)}  dW  (63) 

*0 

and  gives  the  relative  temperature  distribution  through  the  liquid  layer. 

It  Is  evident  that  the  surface  temperature  is  n<y  \ by 

T ■ T + ^ N (A,l } (64) 

It  is  a fairly  comp  Icated  procedure  to  actually  calculate  liquid  layer 
thicknesses  s and  melt  velocities  vffl  as  a function  of  the  laser  absorbed  Irra- 
dlance  aq  and  scale  length  Rp.  In  order  to  carry  out  the  above  computations, 
however,  one  must  be  able  to  calculate  the  vapor  pressure  as  a function  of 
temperature.  The  technique  for  doing  this  Is  given  In  Appendix  A.  Assuming 
this  technique  Is  available,  one  generates  curves  of  s and  vm  as  a function  of 
aq  by: 


1)  Choose  the  coefficient  A and  compute  WA,1).  From  equation  62, 
the  product  (v^s)  Is  thus  set. 

2)  From  equation  Alb  and  the  above  value  of  A compute  the  product 

(soq). 


m 


N (A,  wi 


3)  Equation  64  then  yields  the  surface  temperature  T . 

4)  From  the  results  of  Appendix  A compute  the  reaction  pressure  P . 

5)  From  equation  53c  compute  s. 

6)  Using  the  results  of  steps  1,  2 and  5 compute  aq  and  vffl. 

The  above  procedure  gives  one  point  on  a curve  of  vm  (or  s)  vs. aq  and  the 
above  procedure  Is  repeated  until  a complete  curve  Is  generated. 

For  high  laser  Intensities,  which  lead  to  high  surface  temperatures,  the 
velocity  of  recession  of  the  liquid  surface  due  to  vaporization  vy  Is  not 
negligible  compared  to  vm-  As  noted  previously,  the  vapor  pressure  is  a very 
strong  function  of  the  surface  temperature.  For  this  reason  the  value  of  Rp 
(pressure  scale  length)  Is  generally  less  than  the  laser  spot  size  itself. 

When  large  velocities  vy  occur,  the  liquid-gas  boundary  condition  is  a function 
of  the  radius.  This  causes  all  variables  to  be  dependent  on  r Instead  of 
just  P and  vr.  Section  7.3  considers  the  case  where  significant  mass  is  lost 
In  the  vapor  state. 
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3-  TWO  - DIMENSIONAL  MELT  REMOVAL 


We  next  consider  the  (azimuthal ly  symmetrical)  two-dimensional  melt 
removal  process  taking  Into  account  the  radial  variation  of  the  melt  velo 
city  caused  by  the  radial  variation  of  laser  Intensity.  As  a starting 
point  we  assume  the  melt  velocity  to  have  a parabolic  profile 


m 


mo 


(1 


(65) 


In  general,  higher  order  terms  In  r could  be  included  in  equation  S5.  However, 
for  simplicity,  we  use  only  the  first  two  terms.  In  addition  the  vaporization 
recession  rate  (discussed  In  Appendix  A)  1$  assumed  to  have  tht  form 


v ■ v (1 
v vo  v 


rf) 


(66) 


As  discussed  In  section  7.2,  we  use  a coordinate  system  moving  with  the  local 
melt  boundary  (see  figure  76  ).  Because  of  the  radial  variation  of  vm,  the  melt 
boundary  is  not  strictly  plane. 


Because  of  the  coordinate  system  used  (x  » z - vmt)  the 
tlves  that  occur  In  the  conservation  equations  are  replaced 
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(67) 


With  equation  67,  the  conservation  of  mass  equation  can  be  written 
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Using  the  boundary  conditions  vyi  (r=0) 

r2 

vz  [x*-sl  = V (1  - TO1 


vz(x*0)-!0  and 
vvo  (1  ' 


(68) 


(6S) 


(70) 


(71) 


The  solutions  given  by  equations  70  and  71  are  approximate  since  the 
boundary  condition  vr  (x=o)  =0  is  not  rigorously  satisfied.  However,  as 
can  be  seen  from  the  second  term  of  equation  70,  the  error  is  small. 


i 


Substituting  equations  70  and  71  Into  the  momentum  equation  one  determines 
the  pressure. 


P --  PQ  (1 


(72) 


Again,  as  in  section  7.2,  the  fluid  pressure  is  nearly  constant  through  the 
melt  depth.  Equation  72  has  neglected  terms  of  order  (s/RD)2,  (s/f^)2,  and 
small  terms  Involving  viscosity.  These  neglected  terms  contribute  errors  of 
about  one  percent  if  the  melt  depth  to  radius  ratio  is  less  than  0.1.  This 
assumption  was  satisfied  for  the  results  presented  earlier. 
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The  heat  conduction  equation  can  be  written,  :n  terms  of  the  variable  x 
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Equation  73  has  neglected  terms  in  r and  higher  order  terms. 


The  solution  of  equation  73  Is  written 

|j  » exp  [ -A  B [r.x]]  (yL  + J D [r,x]  exp  [a  B [r,x]]  dx) 

where  the  various  functions  and  coefficients  are  defined  in  Appendix  B. 
Using  those  results  leads  to 

T [r.x]  - D [x]  (tsO  -4>-Tn)  + Tm 

kt 

where  D(0)=0  and  D(-s)=l  satisfy  the  thermal  boundary  conditions. 

At  the  melt-solid  boundary  (x*0)  the  heat  flux  continuity  equation  must 
be  satisfied. 

h ‘ V (|,sns  * ps  Cs  (Tm  • To>  ) 

Equation  76  therefore  shows  that  the  "constant"  yL  has  the  same  radial 
variation  as  does  vm  as  shown  1r  equation  65. 

At  the  liquid-gas  boundary,  at  x»-s,  the  heat  flux  continuity  relation- 
ship leads  to  the  equality  of  thermally  conducted  heat  into  the  liquid  and 
the  net  power  absorbed  q^.  The  quantity  q^  Is  the  amount  of  absorbed  laser 
Irradlance  minus  the  vaporization  energy  carried  away  from  the  surface  and 
given  In  Appendix  A. 
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% ■ »q  - qv 


(77) 
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In  general,  the  scale  length  for  the  variation  of  q is  dependent  on  the 
laser  spot  size  Rg  while  the  scale  size  for  qy,  RyJ  is  usually  significantly 
smaller  than  the  temperature  scale  length  itself. 

r2 

qv = %o  (1  “ r?) 

Rv 

After  a lengthly  series  of  manipulations  one  finds 
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J*'  FJ  [A.*.']  )))  *<* 

The  functions  FQ,  F^ , F2,  Fg  are  again  defined  in  Appendix  B. 

Using  algebra  similar  to  that  used  in  section  7.2  (equating  terms  which 
are  net  a function  of  the  radius)  one  also  determines  the  maximum  (at  the 
center  of  the  laser  spot)  melt,  recession  velocity 

K - qyo)  exp  [-AO  - ^ VJ]  - «\STSF0  [A.V.l] 
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Equation  80  accounts  for  the  fractional  energy  going  into  both  melting 
and  vaporization  and  also  accounts  for  any  radial  heat  flow  in  the  liquid. 
For  very  high  vaporization  rates,  the  liquid  layer  thickness  will  be  small 
(s-*-0)  and  equation  80  reduces  to 


= vmo  (ps  "m  + *s  Cs  <VTo»  W + \0  <81> 

which  is  equivalent  to  the  standard  formulation  of  laser  vaporization  (qVQ 
is  proportional  to  %0Psnv)- 


After  considerable  work  the  liquid  temperature  as  a function  of  radius 
and  depth  is  determined  by  integration  of  equation  74. 
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The  functions  EQ,  Ej,  Eg,  Eg,  GQ,  G-j , Gg,  83,  and  Hj,  Hg,  H3  are  integrals 
dependent  on  the  three  parameters  A,  V and  W.  The  definitions  are  again 
given  In  Appendix  B.  The  maximum  surface  temperature  in  the  center  of  the 
melt  Is  determined  from  equation  82,  for  w » 1,  to  be 

SV 

Tm  + lf  (»s%  + ’s^VV)  Eo  tA.V.U 

T,  - T [o,-*)  5 — (83] 

1 - 4-|-  G [A.V.l] 
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From  equations  75  and  82,  one  also  obtains 


D [x]  = 


svmo  / f #T  T A r m v x-i  , 4 s TsGo  i-A*  v*  " s’] 
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As  described  in  section  7.2,  the  process  of  actually  computing  curves 
showing  the  dependence  of  melt  thickness  s and  melt  velocity  vm0  as  a function 
of  absorbed  laser  irradiance  is  accomplished  by  an  iterative  procedure. 

1)  Select  a numerical  value  for  A as  defined  by  equation  62  which  also 
selects  the  product  s v^. 

2)  Select  a numerical  value  of  the  pressure  scale  length  Rp. 

3)  Select  trial  values  of  s,  Ts,  Ry,  and  V. 

4)  Compute  the  Integrals  G0  and  EQ  which  lead  to  a new  estimate  for  Ts. 

5)  From  the  results  of  Appendix  A compute  the  reaction  pressure  PQ  due 
to  vaporization  of  the  liquid  surface  at  temperature  T$. 

6)  From  Appendix  A,  compute  the  surface  recession  rate  Vy0  due  to 
vaporization. 

7)  From  the  radial  momentum  equation  note  that  the  peak  pressure  PQ 
and  vyQ  are  related  by 


p * RP  vmo  (vmo~vvo^ 
o 2 s* 

8)  From  equation  85  and  step  5 a fourth  order  polynomial  in  s is  obtained 


(85) 


*4  * -r^-  ( (s  VIK>>  (s  \o>  - (s  v»>2)  ■ 0 ,86> 

9)  Equation  86  Is  solved  numerically  for  .*  new  value  of  s by  Newton's 

method.  This  value  of  s and  the  product  s vmo  from  step  I glvns  a 

new  value  of  v_-.  u 

mo  v 

10)  A new  value  of  V * 1 - Is  then  computed. 

vmo 

11)  A new  valye  of  RT  Is  computed  from  an  equation  obtained  by  equating 

1 2 

the  terms  proportional  to  r In  equation  82. 

12)  Steps  4 through  11  are  Iterated  until  all  computed  values  converge. 
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Although  the  above  12  steps  seem  laborious,  they  have  been  carried  out 
for  a range  of  absorbed  laser  Intensities  and  laser  spot  sizes.  The  results 
are  discussed  In  Section  7.1  where  the  melt  thickness,  melt  velocity,  specific 
impulse  I/E,  and  the  specific  mass  removal  efficiency  are  computed  for  aluminum, 
titanium,  and  molybdenum,  and  shown  In  figures  63  through  75. 


4.  TIME  TO  REACH  STEADY  STATE 

The  analyses  given  In  sections  7.2  and  7.3  have  assumed  that  a steady  state 

vaporization  and  melting  situation  extsts.  In  reality*  of  course,  this  steady 

state  takes  a certain  time  to  occur.  If  the  applied  laser  Irradlance  Is  a 

step  function  of  time,  the  temperature  profile  with  a semi-infinite  solid  Is 
(121 

given  by ' 


where  the  surface  temperature  (at  z ■ 0)  Is 


T[..t]-T0**.,  (^y/2  <88> 

It  is  easily  observed  from  equation  87  that  the  heat  penetration  depth 

« is 


For  aluminum,  as  an  example,  the  penetration  depth  is  approximately  0.01  cm 
in  a time  of  100  microseconds. 

Using  the  results  of  equation  88  the  surface  temperature  of  aluminum 

Is  plotted  as  a function  of  time  in  figure  78  assuming  the  Initial  temperature 

was  1OG0K.  The  square-root  dependence  of  the  temperature  on  time  Is,  however, 

valid  only  up  to  the  point  where  significant  melting  and  vaporization  occurs. 

6 2 

For  aq  » 10  watts/cm  this  occurs  In  less  than  20  microseconds.  When  <*q  * 

c p 

103  watts/cm  this  time  Increases  by  a factor  of  100. 

When  the  power  loss  due  to  vaporization  becomes  Important  the  surface 
temperature  T$  Is  a solution  of  the  Integral  equation 
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dx 


In  equation  90,  the  second  term  Is  the  same  as  shown  In  equation  88  while  the 
third  term  shows  the  cooling  effect  caused  by  vaporization.  The  vapor  heat 
loss  qy  Is  described  In  analytical  terms  In  Appendix  A (equation  A23)  and 
depends  explicitly  on  the  vapor  pressure  and  temperature  of  the  surface.  The 
solution  of  equation  90  has  been  solved  numerically  for  aluminum  and  the  results 
are  again  shown  In  figure  78.  The  saturation  in  the  temperature  rise  Is 
clearly  caused  by  vaporization.  The  various  curves  show  how  the  surface 
temperature  falls  after  the  laser  pulse  is  terminated. 


Figure  79  shows  that  the  reaction  pressure  on  the  surface  continues  to 
increase  for  a considerably  longer  time  than  may  be  apparent  from  the  corres- 
ponding slow  increase  of  temperature  shown  in  figure  78.  This  is  indicative 
of  the  very  sensitive  dependence  of  the  vapor  pressure  on  surface  temperature. 
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Figaro  79.  Ruction  Pressure  Due  to  Evaporation  of  Aluminum 
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5.  MELT  REMOVAL  BY  PULSED  LASER  BEAMS 


For  repetitively  pulsed  laser  Interaction  with  targets  It  Is  Important 
to  determine  the  ability  of  a partially  vaporizing  surface  to  "blow"  away 
a thin  liquid  layer  of  molten  metal.  As  demonstrated  in  Sections  7.2  and  7.3 
the  melting  efficiency  o * CW  laser  beams  can  be  relatively  high  If  the 
irradlance  is  not  too  high  (or  the  laser  spot  size  not  too  large).  The 
reaction  pressure  due  to  evaporation  produces  a radial  pressure  gradient  which 
keeps  the  product  of  the  melting  velocity  and  the  melt  layer  thickness  at  a 
relatively  small  value.  Otherwise  the  melting  efficiency  Is  small  and  mass 
removal  by  vaporization  from  the  surface  dominates  over  melt  removal  by 
radial  flow  of  the  liquid.  In  these  large  Intensity  and/or  large  beam  size 
cases  the  melting  efficiency  may  be  Improved  by  superimposing  high  intensity 
pulses  on  a steady,  lower  Intensity  beam.  This  would  keep  the  liquid  layer 
thickness  small  so  that  the  melting  beam  would  be  more  efficient  and  also  at 
the  same  time  supply  pressure  pulses. 

Consider  the  case  of  high  Intensity  laser  beam  pulses  which  cause  a 
pressure  gradient  at  the  target  surface  superimposed  on  a steady  moderate 
intensity  laser  beam  which  melts  the  target  surface  and  determines  the  pulse 
width  and  pulse  repetition  rate  for  maximum  melt  removal.  The  surface 
pressure  as  caused  by  either  LSD  waves  or  vaporization  Is  approximated  by 

P ' P°  ^)2)  (91) 

For  the  part  of  the  flow  which  Is  not  In  the  boundary  layer  near  the 
solid  surface  the  viscosity  terms  In  the  Navler-Stokes  equations  are 
negligible  and  the  momentum  equation  for  the  radial  velocity  component 
becomes 
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(92) 


9vr  ♦ v 

3vr  + v 

3V 

3t  r 

3 r ' z 

az 

1 if. 

pl  ar 


For  the  parabolic  pressure  variation  of  equation  9.1,  the  pressure 
gradient  term  Is 


1 if. 
pl  at* 


r 

2 

T 


In  which  the  characteristic  acceleration  time  Is 


Rp  (^) 


An  Infinite  series  solution  of  this  equation  Is 


(93) 


(94) 


(95) 


This  series  solution  Is  valid  from  * » 0 to  about  * s 1.  Figure  80 
shows  plots  of  the  first  five  and  six  terms  of  the  Infinite  series.  For 
small  times  the  v«av  term  Is  negligible  and  the  radial  velocity  increases 
linearly  with  time.  As  the  velocity  Increases,  the  v*av  term  becomes 
Important  until  at  long  times  a constant  velocity  is  approached  and  the 
time  derivative  term  becomes  negligible. 


V <7 


>-f 


(96) 


A closed  form  expression  which  matches  well  the  solution  from  * * 0 to 
^ ■ 1 and  Is  asymptotic  to  the  constant  long  time  solution  is 

vr'^  (’  + <7)")"^ 


(97) 
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gure  80.  Radial 


with  n > 2.5,  and  Is  shown  as  curve  B in  Figure  80.  A more  approximate 
closed  form  expression  which  matches  the  short  time  and  long  time  solution 
but  gives  slightly  low  velocity  values  for  Intermediate  times 


Vr 

Is  shown  as  curve  C In  Figure  80.  Equation  98  has  the  advantage  that 
It  Is  analytically  Integrable.  It  also  allows  approximately  fqr  the  small 
boundary  layer  which  has  a smaller  average  velocity  thar.'  the,rfest  of  the 
liquid. 

Neglecting  melting  during  the  short  pressure  pulse,  the  time  rate  of 
decrease  of. the  liquid  volume  In  the  cylindrical  melted  region  equals  the 
rate  of  flow  radially  out  of  the  region  through  the  area 


dVol 

ST 


- V * < 


dh 

HF 


(99) 


The  resulting  differential  equation 


has  the  solution  (when  n ■ 2) 


(100) 


2 

h - hQ  exp  [2-2  (l  + ^ ) ] (101) 

T ' 

2 

The  melt  removal  rate  reaches  a maximum  when  t * j t and  then  decreases 
so  that  at  about  t ■ 1.4t  the  rate  approximately  equals  the  rate  at  t - 0.2t. 


For  pulse  durations  which  are  snort  compared  with  the  characteristic 
acceleration  time  x,  the  accelerafor  Is  constant  and  the  solution  for  the 
liquid  depth  Is 

h - h0  exp  [ - (*■)  2 ] (102) 


This  expression  Is  plotted  In  Figure  81.  When  the  pressure  pulse  terminates 
at  t * tj.  the  fluid  continues  to  flow  In  a radial  direction  with  the 

rti 

velocity  vr  = -g-  that  It  had  at  t » tt.  The  differential  equation  for  the 

T 

melt  depth  Is  again 


*2n  r h vr 


2 dh 
*r  Ht 


(103) 


with  constant  v„,  which  has  a solution 
r 


h = hi  exp 


[■*¥] 


7" 


■j  < t < 00 


(104) 


These  expressions  are  also  plotted  in  Figure  81  for  t4  * O.lx,  0.2r  and 
0.3t. 


It  is  seen  that  If  the  pressure  pulse  duration  Is  short  the  melt  doesn't 

get  accelerated  much  and  the  coasting  velocity  is  correspondingly  small. 

Some  combinations  for  effective  melt  removal  are  pressure  pulsing  for  O.lx 

and  then  coasting  for  lOx  or  pressure  pulsing  for  0.3x  and  then  coasting 

for  3x,  etc.  For  Pn  • 1 atmosphere,  Rn  a 1 centimeter  and  liquid  aluminum 

0 3 v 

density  of  2.01  grams  / cm,  the  characteristic  acceleration  time  Is  about 
1 millisecond.  This  relatively  long  time  makes  it  difficult  to  find  a 
realistic  case  where  pressure  pulsing  would  clearly  be  effective. 

Consider  as  some  examples  an  aluminum  target  with  R_  - 1 centimeter 

4 5 « P 

and  steady  absorbed  intensities  of  10  , 10  , and  10  watts  per  square  centi- 
meter. A measure  of  the  times  required  to  reach  steady  state  melting  after 
the  last  pressure  pulse  pushed  away  the  previously  molten  material  can  be 
obtained  by  dividing  the  steady  state  melt  thicknesses  1100pm,  300pm,  and 
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Figure  81.  Melt  Removal  by  Pressure  Pulsing  and  Coasting 
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125ym  by  the  steady  state  melting  velocities  2.3  cm/sec,  14  cm/sec  and 
44  cm/ sec  to  get  48  mi  111  sec,  2.1  millisec  and  0.28  ml  Hi  sec,  respectively. 

As  an  example  let  the  pressure-producing  laser  pulses  have  ten  times  the 
absorbed  intensities  of  the  steady  Intensities  and  f e applied  for  a 
duration  of  C.lx  and  repeated  every  lOx. 

The  pressure  pulses  would  then  have  one  tenth  the  overage  intensity 
cf  the  steady  beam.  Ten  times  the  characteristic  acceleration  times 
would  then  be  22  millisec,  3.8  millisec  and  1.1  Ttillisec  respectively. 

In  only  the  lowest  intensity  case  (1C4  W/cm2  CW  ar.d,  10*  W/cm2  pulsed;  is 
the  coasting  time  required  to  remove  the  less  than  the  melting  time 
to  make  the  melt. 

At  20  pulses  per  second  and  with  a pulse  length  of  200  microseconds 
the  average  Intensity  of  the  pulsed  beam  would  be  4 percent  of  the  steady 
intensity.  In  this  case  the  melting  efficiency  would  only  be  Improved 
slightly  because  the  steady  state  melting  is  already  In  the  efficient  melt- 
ing range.  For  the  middle  intensity  case  (10*  W/cm2  CW,  10*  W/cm2  pulsed), 
pulsing  and  coasting  could  be  used  with  a pulse  duration  of  0.2v  and  a 
coasting  time  of  or  with  an  average  pulsed  intensity  of  40  percent  of  the 
CW  intensity.  For  the  highest  Intensity  case  (10*  W/cm2  CW,  106 7  W/cm2  pulsed) 
pulsing  and  coasting  could  be  used  with  a pulse  duration  of  0.3t  and  a 
coasting  time  of  2.5c  with  an  average  pulsed  Intensity  of  120  percent  of 
the  CW  Intensity.  In  these  latter  cases,  with  relatively  large  proportions 
of  the  total  average  intensity  in  the  pulsed  component,  it  is  unclear 
whether  there  Is  a gain  or  loss  in  melting  efficiency.  The  pulsed  component 
Increases  the  melting  efficiency  of  the  CW  component  but  the  melting  due  to 
the  pulsed  component  Is  less  efficient  than  it  would  be  if  its  average 
intensity  were  applied  steadily.  The  situation  is  similar  for  other  beam 
radii. 

6.  EFFECT  OF  GRAVITY 

Whan  the  force  of  gravity  p 9 Is  superimposed  on  the  radial  pressure 
dp 

gradient  ^ , the  melt  flows  radially  outward  from  a center  which  is 


pl9  Rp 

TFT~ 


ccs  e above  the  center  of  vaporization  pressure  at  r = 0 


in  the  center  of  the  melted  region.  6 is  the  angle  between  the  direction 
normal  to  the  surface  and  the  horizontal  plane.  As  long  as  this  center  of 
flow  Is  Inside  the  melted  region  the  previous  analysis  applies  with  a great- 
er amount  of  the  melt  flowing  out  of  the  crater  or  piling  up  on  the  rim 


on  the  bottom  than  on  the  top.  When  P„  is  less  than  9PL  RP  cos  6,  the 

0 ~i — 

position  where  the  resultant  force  is  zero  is  above  the  melted  region  and 
the  boundary  condition  of  the  flow  going  outward  along  r ■ R is  not 

f 

satisfied  and  the  analysis  Is  modified  and  more  complicated.  This  occurs 

for  PQ  - 0.973  x 10"3Rp  atmosphere/cm  for  aluminum,  PQ  - 1.62  x 10_3Rp 

< -3 

atmosphere/cm  for  titanium  and  PQ  -3.69  x 10  Rp  atmosphere/cm  for 

4 

molybdenum  for  vertical  surfaces.  This  occurs  at  approximately  c*q  * 10 
2 

watts/cm  for  the  various  materials  and  beam  radii. 


SECTION  VIII 


CONCLUSIONS 

The  measurement  of  the  conditions  for  Initiation  of  target  plasmas 
has  been  carried  out  for  both  2.8  and  3.8  microns  wavelength.  Both 
fluence  and  irradlar.ce  (or  temperature  and  irradiance)  thresholds  are 
simultaneously  required  In  order  to  ignite  plasmas.  This  is  evident 
by  comparing  the  theoretical  description  of  these  thresholds  with  the 
wide  range  of  target  materials  and  the  two  wavelengths.  It  Is  evident 
that  the  difference  In  wavelength  (a  factor  of  1.36)  yields  significantly 
different  thresholds  for  our  particular  pulse  shape.  Additional  measure- 
ments and  analysis  In  this  area  would  be  of  great  Importance  in  under- 
standing enhanced  coupling  phenomenology. 

The  technique  of  measuring  the  heat-addition  profile  by  measuring  the 
time  dependence  of  the  temperature  of  the  r « 0 position  has  proven  valuabl 
However,  the  data,  as  analyzed,  are  not  completely  self-consistent.  It 
appears  that  a slngle-Gaussian  description  of  this  heating  profile  Is  not 
adequate.  Because  of  this,  a true  mathematical  inversion  of  the  data 
should  be  carried  out.  Since  this  technique  has  not  been  developed,  either 
electronically  or  mathematically,  additional  work  is  needed.  This  develop- 
ment would  greatly  aid  In  understanding  the  mechanisms  of  enhanced  coupling 
with  either  small  or  large  spots. 

Only  a limited  number  of  measurements  wire  carried  out  showing  the 
dependence  of  mass  loss  of  ablators  with  laser  energy.  These  data  (at 
both  HF  and  OF  wavelength) show  Interesting  limitations  at  both  low  and 
high  delivered  energies.  This  Is  due,  for  low  en?rgies,  to  the  fact  that 
no  vaporization  occurs.  At  high  energies,  on  the  other  hand,  plasma 
production  and  absorption  become  Important. 

One  of  the  goals  of  this  work  was  to  determine  mechanisms  for  melt 
removal  by  Intense  beams.  Experimentally,  the  results  were  not  Impressive 
since  molten  metals  were  not  directly  removed  In  large  amounts  from  only 
partially  melted  surfaces.  The  surface  pressures  do,  however,  push  the 


molten  material  through  a completely  melted  surface.  This  area  of 
Investigation  obviously  requires  more  experiments  over  larger  beam  spot 
sizes  and  correspondingly  larger  laser  energies.  (The  above  experiments 
were  done  with  25  joules  of  C02  laser  energy).  Theoretical  analysis  of 
melt  removal  by  a partially  vaporizing  surface  leads  to  the  expectation 
that  high  specific  mass  removal  efficiencies  can  be  obtained,  even  In  the 
absence  of  an  external  wind  or  gravity.  This  results  from  the  gradient 
of  the  vapor  pressure  causing  a radial  flow  of  the  liquid  metal  surfaces. 
The  analysis  developed  should  be  improved,  however,  and  additional  data 
presented. 
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APPENDIX  A 

VAPOR  PRESSURES  AND  EVAPORATION  RATES 

When  two  different  phases  (solid,  liquid  or  gas)  of  a substance  are 

In  equilibrium,  the  temperature  T,  pressure  P and  chemical  potential  p 
/?4) 

must  be  equal  v 

{P*T>  * n2  (P.T)  (Al) 

so  that  the  pressure  and  temperature  of  phases  In  equilibrium  can  be 
expressed  as  functions  of  each  other.  When  one  of  these  quantities  (P.T) 

Is  given  the  other  Is  completely  determined.  Differentiating  the  chemical 
potential  equation  with  respect  to  the  temperature  yields 


3yl  ^1  dP  3y2  8y2  dP 

rr+  ar  ar’  ft  + wxr 


(A2) 


Then,  using  the  Identity 

dp  3 - s dT  + v dP  (A3) 

In  which  s Is  the  molecular  entropy  and  v Is  the  molecular  volume,  and  the 
expression  for  the  latent  heat  0 of  the  transition 

fl  ■ T(Sg  - Sj)  (A4) 

gives  the  Clausius  - Clapeyron  equation 

dP  _ n 

ar  nv2"-vp  (as) 

This  aquation  gives  the  change  of  pressure  with  change  of  temperature  along 
the  phase  equilibrium  curve. 

*4.  L4nd«u,  t.  D.  and  vtfstoix,  t.  Statistical  Physics.  Addison -Wesley 
Publishing  Company,  Inc.  1958. 

m 


The  molecular  volume  of  a gas  v2  ■ p~  Is  generally  much  larger  than 

the  molecular  volume  of  the  condensed  phase  so  for  the  case  of  a solid  or 
liquid  In  equilibrium  with  Its  vapor,  v^  can  be  neglected  and  the  Clausius- 
Clapeyron  equation  becomes 

dP  „ nP 

ar  jjs  (ac 

If  the  latent  heat  of  vaporization  were  constant  the  vapor  pressure  would 
vary  exponentially  with  the  reciprocal  of  the  temperature. 


P ■ PQ  exp  {- 


Figure  A1  shows  the  specific  heats  of  aluminum,  titanium  and  molybdenum. 
The  specific  heats  of  the  liquids  are  constant.  The  specific  heat  of 
aluminum  gas  near  the  boiling  temperature  Tb  s 2766. 8°K  (at  which  the  vapor 
pressure  Is  one  atmosphere)  Is  very  nearly  constant  so  that  the  latent  heat 
of  vaporization  Is  a linear  function  of  the  temperature. 

Thus 

^ " ^e02  " e01^  + ^p2  * T (A* 


where 


e02  ’ e01 


- (Cp2 


C,)  T„ 


Is  the  latent  heat  of  vaporization  extrapolated  to  zero  temperature.  With 
equation  M the  solution  of  equation  A6  becomes 


(0 


- C,) 


^02^01  ^ 


^ Y - y ^ J Atmospheres  (A9) 


TEMPERATURE  (1,000  K) 

Figure  AT.  Specific  Heats 


Using  the  .1ANAF  tables,  reference  25,  one  obtains  (for  aluminum) 


1.32 


(A10) 


38,620°  K 


The  vapor  pressure  of  aluminum,  found  from  equations  A9  and  A10,  Is  plotted 
In  Figure  A2.  The  x symbols  are  from  Reference  2b,  the  + symbols  are  from 
reference  26  and  the  other  symbols  are  from  other  references. 


From  equation  A9,  one  finds  at  T a T^  * 2767°K  that 


dP 

3T 


12.64  j . 


(All) 


Thus,  from  the  assumption  of  Section  7.2  for  the  vapor  pressure  variation 
Kith  transverse  distance  [ P * PQ  (1 -r2/P^  )]  and  assuming  that 

T - T$  (l-r2/R?)  (A12) 


one  finds  that  (for  aluminum) 
kp  - 0.28  Rt 


(A12) 


At  other  temperatures,  one  finds  from  equations  A5  and  A8, 


I Aluminum 


( 38620 

V nr 


(A13) 


25.  JANAF  Themochemi cal  Data,  The  Dow  Chemical  Company,  Thermal  Laboratory, 
midland,  Michigan,  1968. 

26.  Handbook  of  Chemistry  and  Physics,  p.  0-142  of  49th  Edition,  The 
chemical  Rubber  to.,  Cleveland,  Ohio,  1969. 


171 


As  expected  the  scale  length  of  the  pressure  profile  Is  considerably 
smaller  than  the  temperature  scale  length. 

The  solution  of  the  Clauslus-Clapuyron  equation  for  titanium  Is 
somewhat  more  complicated  than  for  aluminum  due  to  the  Increase  of 
specific  heat  of  the  gas  with  temperature.  As  seen  from  Figure  A!  one 
may  approximately  write 

Cp2  -VCP  <T-V 


a 


e02  " e01 


+ (cp  ‘C1  -cp  Vr  + cp  7 


(A14) 


e02  “ e01 


AHU  - (Cn  -C,)Tk  + -£ 


T2 


v '-p  w-j ' 'b  T T 


Using  equations  A14  and  A6,  and  numerical  parameters  found  In  reference  23, 
one  determines  (for  titanium) 


P " (355T-)  exp[  ”56280  (y  - ‘ ] Atmospheres 

Equation  A15  Is  plotted  In  Figure  A2. 

Again  one  can  determine  (by  the  same  technique  as  used  for  aluminum) 
the  ratio  of  temperature  and  vapor  pressure  scale  lengths  for  a titanium 
target  material. 


(A15) 


rt 


„ ( 56280  ? co  + T \ 

v-t-  " 2 68  + im) 
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(A16; 


Titanium 


For  titanium  then*  one  finds  for  Ts  - 3591°K  and  P ■ lAtm  that  Rj/Rp  c 3.77. 

Raising  the  titanium  surface  temperature  to  5000°K  (P  ■ 55.7  Atm)  leads  to 
Ry/Rp  - 3.21. 
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For  molybdenum,  the  same  type  analysis  as  used  for  titanium  applies 
as  can  be  seen  from  Figure  Al.  The  vapor  pressure  Is  found  to  be 

p ‘(tot?)  6 36  •x|)  [-  83590  (t  - tot?)  + Irar  ] 

' ’ (A17) 

Equation  A17  Is  again  plotted  In  figure  A2.  The  scale  lengths  for  temperature 
and  pressure  are  again  found  to  be  related  by 


83690 


6.36  + 


T35S 


(A18) 


At  atmospheric  pressure,  T * 4919  K,  one  obtains  fU  * 3.8  Rp. 


When  calculating  forces  on  laser  vaporized  surfaces,  it  is  worthwhile 
to  recall  several  elementary  principles.  First  of  all,  the  reaction  pressure 
due  to  vaporization  Is  one  half  the  corresponding  vapor  pressure  as  determined 
above.  This  occurs  because  of  the  momentum  reversal  of  molecules  incident 
on  a noncondensing  surface.  The  rate  of  molecular  transfer  n through  the 
vaporizing  surface  Is  given  by 

• / «/■►  \V2  o 


• / kt  V/c  2 

n * nQ  ( ) particles/cm  /sec 


(A19) 


where  nQ  Is  the  gas  density  and  m is  the  mass  of  the  vapor  molecule.  Since 
Py  ■ nkT,  equation  A19  can  be  rewritten 


(A20) 


(2mkt) 

Since  the  equilibrium  rate  of  transfer  of  particles  across  a surface  in  one 
direction  is  the  same  as  In  vacuum,  equation  A20  leads  to  the  recession  rate 
vv  of  the  liquid  surface.  From  continuity  of  mass 


"a 


P / X 1/2 

■ M-3  sf  (t;-)  a"/** 


(A21 ) 
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In  equation  A21,  A Is  the  atomic  weight,  NQ  Is  Avogadro’s  number,  and 
Is  the  mass  density  of  the  liquid. 


The  kinetic  energy  transfer  through  the  liquid-gas  surface  q^  Is  similarly 
found  to  be 


• ? 
qk  * 2kT  n watts/cm 


(A22) 


The  total  power  loss  from  the  surface  qy  Includes  the  heat  of  vaporization 
so  that 


qy  « n(2kT  + U) 

PV(2T,  + £)  , 

* 369.  - — watts/cm 

( ATS  )1/2 

where  Py  Is  given  In  atmospheres. 


(A23) 


From  Reference  25  we  obtain 


a 

k « 38,620  - 1.32  T 

A1  s 


ft 

F 


- 56,280 


2.68  Ts 


(A24) 


- 83,590  - 6.36  T$ 


+ 


vm 


Thus  we  see  that  the  heat  of  vaporization  terms  are  considerably  larger  than 
the  corresponding  kinetic  energy  terms  of  equation  A23. 


APPENDIX  B 


FUNCTIONS  RELATED  TO  THE  TWO-DIMENSIONAL 
HEAT  FLOW  SOLUTION 

Section  7.3  describes  the  two-dimensional  temperature  distribution  in 
terms  of  several  functions  and  parameters. 
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Including  only  the  first  two  terms  in  the  expansion  of  equation  B4  yield? 


■M 
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The  functions  FQ,  F1 , Fg,  and  F3  are  numerically  determined  and  are 


defined  by 


w 


F„  [A,V,w]  = Jw" 
0 


dw 


(B6) 


The  integrations  of  equation  B6  were  carried  out  using  the  trapezoidal 
3 

rule  and  10  steps. 


The  function  En>  Gn  and  Hn  (where  n = 0,  1,  2,  3)  which  define  the 
temperature  profile  (see  equation  82)  are  defined  by 

w 

En  [ A*  V’  W]  * jwn  exp  [a(w  - j Vw3)]  dw  (B7) 

o 

w 

Jwn  Fq  [A,  V,  w]  exp  [a(w  - 4 Vw3)  ] dw  (B8) 
o 

w 

jFn  [A,  V,  w]  exp  [a(w  - 1 Vw3)]  ciw  (B9) 

o 

In  equations  B8  and  R9  the  functions  Fn  are  defined  by  equation  B6. 


and 

H„[A,  V,  w]  - 


G„[A,  V,  M]  = 
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